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Abstract 

Aspects  of  the  circulation  and  convection  in  the  Irminger  Sea  are  investigated  using  a  vari¬ 
ety  of  in-situ,  satellite,  and  atmospheric  reanalysis  products. 

Westerly  Greenland  tip  jet  events  are  intense,  small-scale  wind  phenomena  located  east 
of  Cape  Farewell,  and  are  important  to  circulation  and  convection  in  the  Irminger  Sea. 
A  climatology  of  such  events  was  used  to  investigate  their  evolution  and  mechanism  of 
generation.  The  air  parcels  constituting  the  tip  jet  are  shown  to  have  a  continental  origin, 
and  to  exhibit  a  characteristic  deflection  and  acceleration  around  southern  Greenland.  The 
events  are  almost  invariably  accompanied  both  by  a  notable  coherence  of  the  lower-level 
tip  jet  with  an  overlying  upper- level  jet  stream,  and  by  a  surface  cyclone  located  in  the  lee 
(east)  of  Greenland.  It  is  argued  that  the  tip  jet  arises  from  the  interplay  of  the  synoptic- 
scale  flow  evolution  and  the  perturbing  effects  of  Greenland’s  topography  upon  the  flow. 

The  Irminger  Gyre  is  a  narrow,  cyclonic  recirculation  confined  to  the  southwest  Irminger 
Sea.  While  the  gyre’s  existence  has  been  previously  documented,  relatively  little  is  known 
about  its  specific  features  or  variability.  The  mean  strength  of  the  gyre’s  circulation  be¬ 
tween  1991  and  2007  was  6.8  ±1.8  Sv.  It  intensified  at  a  rate  of  4.3  Sv  per  decade  over  the 
observed  period  despite  declining  atmospheric  forcing.  Examination  of  the  temporal  evo¬ 
lution  of  the  LSW  layer  thickness  across  the  Irminger  Basin  suggests  that  local  convection 
formed  LSW  during  the  early  1990s  within  the  Irminger  Gyre.  In  contrast,  LSW  appeared 
outside  of  the  gyre  in  the  eastern  part  of  the  Irminger  Sea  with  a  time  lag  of  2-3  years, 
consistent  with  transit  from  a  remote  source  in  the  Labrador  Sea. 

In  the  winter  of  2007-08  deep  convection  returned  to  both  the  Labrador  and  Irminger 
seas  following  years  of  shallow  overturning.  The  transition  to  a  convective  state  took  place 
abruptly,  without  going  through  a  preconditioning  phase,  which  is  contrary  to  general  ex¬ 
pectations.  Changes  in  the  hemispheric  air  temperature,  tracks  of  storms,  flux  of  freshwater 
to  the  Labrador  Sea,  and  distribution  of  pack  ice  all  conspired  to  enhance  the  air-sea  heat 
flux,  resulting  in  the  deep  overturning. 

Thesis  Supervisor:  Robert  S.  Pickart 

Title:  Senior  Scientist,  Woods  Hole  Oceanographic  Institution 
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Chapter  1 


Introduction 


The  Irminger  and  Labrador  seas  form  the  two  deep  basins  in  the  western  subpolar  North 
Atlantic  (Fig.  1-1).  Wintertime  convection  in  this  region  leads  to  the  formation  of  an  in¬ 
termediate  water  mass  called  Labrador  Sea  Water  (LSW),  which  ventilates  a  vast  area 
of  the  North  Atlantic  and  impacts  the  stratification  of  the  water  column.  The  process  of 
convection  has  important  climatic  impacts:  Heat  is  released  from  the  ocean  to  the  atmo¬ 
sphere,  which  energizes  the  North  Atlantic  storm  track  and  contributes  to  the  net  poleward 
heat  transport  by  the  ocean-atmosphere  system.  CO2  dissolved  near  the  surface  is  also 
sequestered  at  depth  by  this  process. 

1.1  Circulation  and  ventilation  of  the  subpolar  gyre 

From  its  formation  region  in  the  northwestern  subpolar  gyre,  LSW  spreads  out  and  influ¬ 
ences  a  large  area  of  the  North  Atlantic  and  beyond  (e.g.  Talley  and  McCartney,  1982;  Stra- 
neo  et  al,  2003;  Yashayaev  et  al,  2007a).  Internal  cyclonic  recirculations  in  the  Labrador 
and  Irminger  seas  maintain  an  exchange  of  LSW  between  these  two  basins.  Along  interior 
pathways  LSW  spreads  to  the  eastern  North  Atlantic,  and  LSW  is  also  a  direct  component 
of  the  Deep  Western  Boundary  Current  (DWBC),  by  which  it  is  exported  equatorward  and 
out  of  the  subpolar  North  Atlantic.^  The  DWBC  is  the  deep  limb  of  the  Atlantic  Merid- 

^The  recent  study  of  Bower  et  al.  (2009)  suggests  that  pathways  in  the  interior  basin  may  be  as  important 
as  the  DWBC  for  the  equatorward  export  of  LSW  south  of  the  Grand  Banks. 
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Figure  1-1:  Topography  of  the  subpolar  North  Atlantic  from  the  Etopo2  2-minute 
base. 
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ional  Overturning  Circulation  (MOC),  which  is  of  critical  importance  to  the  climate  of 
northwestern  Europe  (e.g.  Vellinga  and  Wood,  2002).  The  other  components  of  the  DWBC 
besides  LSW  are  the  overflow  waters  from  the  Nordic  seas:  Denmark  Strait  Overflow  Water 
(DSOW)  and  Iceland-Scotland  Overflow  Water  (ISOW).  As  the  overflow  waters  descend 
from  the  sills  of  the  Greenland-Iceland-Scotland  Ridge,  ambient  water  is  entrained.  Pri¬ 
mary  among  the  entrained  waters  is  LSW  (Fogelqvist  et  al,  2003),  and  in  this  capacity 
LSW  participates  also  in  the  deepest  part  of  the  DWBC.  The  three  water  masses  that  con¬ 
stitute  the  DWBC  (LSW,  ISOW,  and  DSOW)  are  collectively  referred  to  as  North  Atlantic 
Deep  Water  (NADW).  The  properties  of  this  composite  water  mass  are  set  in  the  subpolar 
North  Atlantic,  and  it  is  found  throughout  the  World  Ocean  (Johnson,  2008). 

The  North  Atlantic  subpolar  gyre  circulates  cyclonically  between  about  50°  and  65°N. 
According  to  Sverdrup  theory,  the  gyre  consists  of  a  mass-balanced  slow,  northward  inte¬ 
rior  flow  and  a  swift,  southward  western  boundary  current  driven  by  positive  wind  stress 
curl.  In  practice,  the  subpolar  gyre  circulation  appears  to  be  more  complicated,  and  does 
not  obey  simple  Sverdrup  dynamics  (Bryan  et  al,  1995;  Treguier  et  al,  2005).  A  schematic 
of  the  boundary  currents  in  the  Irminger  Sea  (Fig.  1-2)  shows  warm,  salty  Atlantic-origin 
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water  flowing  northward  along  the  Reykjanes  Ridge  in  the  Irminger  Current  (IC).  After 
a  bifurcation  south  of  the  Denmark  Strait,  where  a  smaller  branch  continues  northward 
through  the  strait  to  form  the  Icelandic  Irminger  Current,  the  bulk  of  the  current  recircu¬ 
lates  to  the  south  to  flow  alongside  the  East  Greenland  Current  (EGC),  which  carries  cold, 
fresh  Arctic-origin  water  equatorward  near  the  shelfbreak.  On  the  inner  shelf  the  East 
Greenland  Coastal  Current  (EGCC),  predominantly  a  bifurcated  branch  of  the  EGC,  ad- 
vects  a  combination  of  Arctic-origin  water  and  coastal  runoff  to  the  south  (Sutherland  and 
Pickart,  2008).  Seaward  of  the  shelfbreak,  on  the  mid  to  deep  continental  slope,  resides  the 
DWBC. 


30W  25W_ 

I  ,'■]  •  A  j.  4  rt' 


Eigure  1-2:  The  boundary  currents  of  the  Irminger  Sea,  from  Pickart  et  al.  (2005).  The  acronyms 
are:  IC  =  Irminger  Current;  EGC  =  East  Greenland  Current;  EGCC  =  East  Greenland  Coastal  Cur¬ 
rent;  DWBC  =  Deep  Western  Boundary  Current.  Red  (blue)  arrows  indicate  upper-layer  transport 
of  warm  (cold)  water,  and  the  numbers  are  transport  estimates  in  Sverdrups  (1  Sv  =  10^  rn^/s).  The 
500,  1000,  2000  and  3000  m  isobaths  are  plotted. 


1.2  Convection  in  the  Irminger  Sea 

Contrary  to  the  implications  of  its  name,  a  growing  body  of  evidence  suggests  that  ESW  can 
be  formed  in  the  Irminger  Sea  as  well  as  in  the  Labrador  Sea  (e.g.  Pickart  et  al,  2003a,b; 
Bacon  et  al,  2003;  Ealina  et  al.,  2007;  Vage  et  al,  2008).  This  idea  is  not  new.  The  no- 
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tion  of  deep  convection  in  the  Irminger  Sea  was  proposed  almost  100  years  ago  by  Nansen 
(1912),  who  presented  a  limited  number  of  hydrographic  stations  with  matching  surface 
and  deep  properties.  Despite  further  evidence  of  deep  convection  in  the  Irminger  Sea  (e.g. 
Wattenberg,  1938),  this  idea  was  largely  forgotten  once  deep  convection  was  directly  ob¬ 
served  in  the  Labrador  Sea.^  The  idea  was  recently  revisited  by  Pickart  et  al.  (2003b), 
who  showed  that  the  observed  lateral  distribution  of  LSW  in  the  western  subpolar  North 
Atlantic  was  inconsistent  with  a  single  source  in  the  Labrador  Sea.  Pickart  et  al.  (2003b) 
pointed  out  that  conditions  favorable  for  open  ocean  convection  existed  in  the  southwestern 
Irminger  Sea  as  well.  In  particular,  during  a  string  of  severe  winters  in  the  early  1990s  the 
atmospheric  forcing,  together  with  the  oceanic  preconditioning  and  cyclonic  circulation, 
satisfied  the  criteria  for  deep  convection  as  outlined  by  Marshall  and  Schott  (1999). 

In  contrast  to  the  Labrador  Sea,  the  conditions  leading  to  deep  convection  in  the  Irminger 
Sea  are  to  a  greater  extent  governed  by  small-scale  atmospheric  and  oceanic  phenomena. 
Consider  first  the  atmospheric  conditions.  Due  to  the  impact  of  the  high  Greenland  topog¬ 
raphy  on  low  pressure  systems  that  pass  nearby,  causing  a  variety  of  intense,  small-scale 
wind  phenomena,  the  region  near  Cape  Farewell  (the  southern  tip  of  Greenland)  is  the 
windiest  location  of  the  World  Ocean  (Sampe  and  Xie,  2007).  Of  particular  importance 
to  wintertime  convection  in  the  southwestern  Irminger  Sea  is  the  westerly  Greenland  tip 
jet  (see  Figure  1-3,  Pickart  et  al,  2003a;  Vage  et  al,  2008).  Tip  jets  are  narrow,  intermit¬ 
tent  wind  patterns  east  of  Cape  Farewell,  whereby  strong  winds  help  elevate  wintertime 
air-sea  heat  fluxes  to  levels  comparable  to  those  in  the  Labrador  Sea.  Although  tip  jet 
events  are  short-lived  (up  to  3  days,  Pickart  et  al,  2003a),  they  occur  repeatedly  through¬ 
out  the  winter  months.  For  example,  during  the  winter  of  1993-94  Pickart  et  al  (2003a) 
documented  10  westerly  tip  jet  events. 

The  small-scale  oceanic  feature  that  is  important  to  deep  convection  in  the  Irminger  Sea 
is  a  closed,  cyclonic  recirculation  gyre  in  the  western  part  of  the  basin  called  the  Irminger 
Gyre  (Lavender  etal,  2000) .  Water  parcels  inside  the  gyre  are  trapped  and  isolated,  thereby 
subject  to  the  strong  atmospheric  forcing  for  a  longer  period  of  time.  The  gyre  circulation 

^The  historical  perspective  of  the  oceanographic  community’s  evolving  view  on  convection  in  the  subpolar 
North  Atlantic  is  detailed  in  Pickart  et  al.  (2008) 
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Wind  speed  (m/s) 


Figure  1-3:  Satellite  seatterometer  (QuikSCAT)  wind  speed  (m/s,  eolor)  and  vectors  showing  a  tip 
jet  on  5  December,  2002. 

also  promotes  retention  of  previous  winters’  convected  products,  which  acts  to  reduce  the 
stratification  of  the  water  column.  Finally,  isopycnal  doming  inside  the  gyre  brings  the 
weakly  stratified  water  closer  to  the  surface.  In  an  integrated  sense  a  series  of  westerly  tip 
jet  events  over  the  course  of  a  winter  have  a  strong  effect  on  the  evolution  of  the  mixed 
layer  of  the  isolated  water  column  near  the  center  of  the  gyre  (Vage  et  al,  2008). 

1.3  Effect  of  atmospheric  variability  on  the  circulation 

The  subpolar  gyre  is  strongly  affected  by  interannual  atmospheric  variability.  The  dominant 
mode  of  variability  in  this  region  is  the  North  Atlantic  Oscillation  (NAO,  e.g.  Hurrell, 
1995).  The  NAO  is  described  by  an  index  computed  from  the  normalized  sea  level  pressure 
difference  between  the  regions  of  the  Icelandic  Low  and  the  Azores  High.  Winters  of  a 
high  NAO  index  are  typically  characterized  by  stronger  westerly  winds,  enhanced  air-sea 
buoyancy  fluxes,  and  deeper  convection.  Modulated  by  significant  year  to  year  variability, 
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Figure  1-4:  Topography  of  southeast  Greenland  and  the  Irminger  Sea  from  the  Etopo2  2-minute 
elevation  data  base.  The  gray  lines  are  contours  of  objectively  mapped  geostrophic  pressure  at  700 
db,  from  Lavender  et  al.  (2000).  The  closed  contours  in  the  western  Irminger  Sea  reveal  the  location 
of  the  Irminger  Gyre. 

the  NAO  index  evolved  from  strongly  negative  values  in  the  mid-1960s  to  strongly  positive 
values  in  the  mid-1990s.  Towards  the  end  of  this  period  the  eoldest  and  densest  LSW 
on  reeord  was  formed  in  large  amounts  (Yashayaev,  2007).  This  trend  abruptly  halted 
in  the  winter  of  1995-96,  and  the  deereased  levels  of  wind  and  air-sea  buoyaney  foreing 
during  the  majority  of  the  subsequent  winters  have  led  to  a  signifieant  reduetion  of  LSW 
formation  (e.g.  Lazier  et  al,  2002;  Yashayaev,  2007). 

The  effeet  of  the  varying  NAO  and  LSW  formation  rates  on  the  subpolar  gyre  was 
nieely  illustrated  by  Hakkinen  and  Rhines  (2004).  Their  analysis  of  sea  surfaee  height 
anomalies  in  the  North  Atlantie  revealed  a  spin-down  of  the  subpolar  gyre  in  eoneert  with 
the  persistent  deeline  of  the  NAO  index  following  the  extreme  positive  state  in  the  mid- 
1990s.  Considering  hydrographie  data  and  the  deelining  buoyaney  foreing,  they  linked  the 
weakening  gyre  eireulation  to  redueed  levels  of  eonvective  activity.  During  this  period  the 
gyre  also  contracted  laterally  (Hatun  et  al,  2005;  Bersch  et  al,  2007).  Similarly,  Curry 
and  McCartney  (2003)  found  that  the  gyre  circulation  reflected  a  time  integral  of  the  NAO 
index,  to  a  large  extent  governed  by  dense  water  formation  in  the  subpolar  North  Atlantic. 
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1.4  Convective  activity  and  the  MOC 


Coupled  ocean-atmosphere  models  integrated  over  long  periods  of  time  often  display  a  co¬ 
herence  between  LSW  formation  and  MOC  intensity^,  with  a  reduced  convective  activity 
in  the  subpolar  North  Atlantic  linked  to  a  decline  in  the  MOC  (e.g.  Jungclaus  et  al,  2005; 
Stouffer  et  al,  2006).  The  performance  of  such  coupled  climate  models  used  in  the  In¬ 
tergovernmental  Panel  on  Climate  Change  (IPCC)  Fourth  Assessment  Report  (AR4)  was 
evaluated  in  the  Labrador  and  Irminger  seas  by  de  Jong  et  al.  (2009).  They  found  that  in 
general  the  simulations  were  less  than  satisfactory,  with  model  biases  exceeding  the  ob¬ 
served  range  of  variability  by  orders  of  magnitude.  The  overall  tendencies  of  the  coupled 
models  were  to  have  a  low  surface  salinity  and  high  intermediate  and  deep  salinities,  which 
led  to  increased  stratification  in  the  upper  1000  db  of  the  models.  Cold  surface  water  above 
warm  intermediate  and  deep  layers  in  most  models  had  a  negative  influence  on  the  stability, 
but  far  from  sufficient  to  compensate  for  the  stabilizing  effect  of  the  salinity.  As  a  result  the 
convective  activity  was  confined  to  a  thin  upper  layer.  Another  group  of  models  had  high 
salinity  and  potential  density  biases  over  the  entire  water  column,  which  led  to  excessively 
deep  convection.  Several  of  the  models  exhibited  stronger  convection  in  the  Irminger  Sea 
than  in  the  Labrador  Sea,  and  formation  of  LSW  was  in  general  poorly  represented  in  the 
coupled  models.  Model  drift  often  became  a  considerable  problem,  in  particular  for  models 
that  required  a  long  spin-up  time.  Since  the  MOC  and  other  aspects  of  climate  variability 
are  highly  sensitive  to  the  hydrographic  conditions  and  circulation  of  the  western  subpolar 
North  Atlantic,  these  low-resolution  simulations  must  be  interpreted  with  caution.  While 
high  resolution  is  required  for  a  realistic  simulation  of  the  circulation  field,  a  redeeming 
feature  of  the  coupled  models  is  that  variability  consistent  with  observed  conditions  can  be 
simulated  with  lower  resolution  (Boning  et  al,  2006). 

A  number  of  high-resolution  hindcast  simulations  of  the  subpolar  North  Atlantic  have 
been  performed  over  the  last  few  years.  While  more  realistic  than  the  low-resolution  sim¬ 
ulations,  such  regional  studies  of  the  subpolar  gyre  circulation  using  general  circulation 

^The  overturning  circulation  may  be  considered  in  both  depth  space  and  density  space,  involving  sinking 
and  water  mass  transformation,  respectively  (e.g.  Pickart  and  Spall,  2007).  Unless  otherwise  noted,  results 
are  discussed  in  depth  space. 
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models  (GCMs)  are  also  challenging  for  a  number  of  reasons,  including  the  importance  of 
buoyancy  forcing  on  the  dynamics  of  the  system,  the  low  stratification,  and  the  difficulty 
of  maintaining  the  small  temperature  and  salinity  contrasts  between  the  water  masses  (e.g. 
Treguier  et  al,  2005;  Maine  et  al,  2008).  It  was  recently  suggested  that  with  a  current 
state-of-the-art  global  GCM  it  is  not  possible  to  achieve  both  dynamic  and  thermodynamic 
fidelity,  in  other  words  realistically  simulating  both  the  intensity  of  the  circulation  and  the 
water  mass  characteristics  at  the  same  time  is  very  difficult  (Yaeger  and  Jochum,  2009). 
Though  directly  applicable  only  to  their  model,  it  may  have  a  bearing  also  on  other  simu¬ 
lations. 

A  shared  feature  also  among  many  of  these  high-resolution  GCMs  is  a  link  between 
convective  activity  in  the  subpolar  gyre  and  the  Atlantic  MOC  (e.g.  Eden  and  Willebrand, 
2001;  Boning  etal,  2006;  Deshayes  and  Frankignoul,  2008).  Closer  inspection  of  some  of 
these  models  reveal  that  even  though  they  tend  to  simulate  the  gross  features  of  the  subpo¬ 
lar  gyre  circulation  reasonably  realistically,  important  features  of  the  observed  conditions 
such  as  the  magnitude  of  the  overflows  and  the  location  and  intensity  of  convection  are 
problematic.  For  example,  in  the  20  km  resolution  model  of  Deshayes  and  Frankignoul 
(2008)  convection  in  the  Irminger  Sea  was  too  intense  and  occurred  in  the  northeastern  part 
of  the  basin  (Deshayes  et  al,  2007),  while  it  is  well  established  that  the  deepest  convection 
occurs  in  the  southwestern  part  of  the  basin  (Pickart  et  al,  2003b;  Centurioni  and  Gould, 
2004).  The  location  matters.  The  close  proximity  of  the  model’s  Irminger  Sea  convec¬ 
tive  region  to  the  boundary  current  led  to  rapid  export  of  newly  convected  water  out  of  the 
basin  via  the  DWBC,  which  resulted  in  an  excessive  influence  on  the  circulation.  Despite 
these  shortcomings,  Deshayes  and  Frankignoul  (2008)  argue  that  the  Atlantic  MOC  can 
to  first  order  be  considered  an  integrator  of  NAO-modulated  atmospheric  forcing,  with  a 
particularly  strong  link  to  convection  in  the  Irminger  Sea. 

Data  assimilation,  the  combination  of  observations  and  model  output,  is  a  common 
practice  in  the  atmospheric  sciences.  This  approach  has  recently  been  applied  also  in  phys¬ 
ical  oceanography  (an  example  of  which  is  the  consortium  Estimating  the  Circulation  and 
Climate  of  the  Ocean  (ECCO),  Wunsch  and  Heimbach,  2007),  where  it  is  referred  to  as 
“state  estimation”.  Essentially  the  state  estimate  is  computed  from  a  GCM  whose  initial 
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and  boundary  conditions  are  subject  to  a  least  squares  adjustment  in  order  to  bring  the  re¬ 
sult  into  agreement  with  observations  within  estimated  errors.  State  estimation  provides  a 
promising  avenue  to  synthesize  diverse  ocean  observations  and  estimate  quantities  that  are 
difficult  to  infer  from  observations  alone,  such  as  the  MOC. 

The  evaluation  of  state  estimates’  realism  and  usefulness  compared  to  GCMs  is  only  be¬ 
ginning  (Schott  and  Brandt,  2007).  Presently  available  computational  resources  dictate  that 
the  resolution  remains  relatively  low  (currently  1  degree  for  the  ECCO  project),  and  this 
is  likely  to  be  a  problem  for  the  foreseeable  future  (Wunsch  and  Heimbach,  2007).  Com¬ 
parison  between  assimilated  and  directly  observed  hydrographic  properties  in  the  Irminger 
Sea  is  not  very  convincing  (de  Jong  et  al,  2009).  Relative  to  observations,  the  mean  ECCO 
vertical  profile  was  significantly  offset  in  terms  of  both  temperature  and  salinity,  and  identi¬ 
fication  of  distinct  water  masses  was  not  possible.  Excessive  convection  was  a  contributing 
factor;  the  water  column  was  regularly  mixed  all  the  way  to  the  bottom  in  both  the  Eabrador 
and  the  Irminger  seas  (Schott  and  Brandt,  2007;  de  Jong  et  al,  2009). 

1.5  Climatic  impacts  of  LSW  formation 

The  importance  of  dense  water  formation  in  the  northern  North  Atlantic  is  well  illustrated 
by  the  meridional  heat  transport  in  the  Atlantic,  which  is  northward  at  all  latitudes,  includ¬ 
ing  south  of  the  equator.  The  magnitude  of  the  heat  transport  has  been  estimated  from  a 
number  of  hydrographic  transects  and  inverse  calculations  in  the  vicinity  of  25°N  to  be 
about  1.3  PW  (1  PW  =  10^^  W,  e.g.  Hall  and  Bryden,  1982;  Ganachaud  and  Wunsch,  2000; 
Talley,  2003).  At  this  latitude  the  oceanic  heat  transport  is  comparable  in  magnitude  to 
the  atmospheric  heat  transport,  while  outside  of  the  tropics  the  atmosphere  dominates  the 
meridional  heat  transport  (Trenberth  and  Caron,  2001).  The  Atlantic  heat  transport  may  be 
divided  into  three  separate  components;  shallow,  intermediate,  and  deep  overturning,  where 
the  latter  two  are  associated  with  ESW  and  NADW  formation,  respectively.  Prom  such  an 
analysis  in  density  space  Talley  (2003)  determined  that  the  three  components  were  of  sim¬ 
ilar  magnitude.  On  the  other  hand,  a  repeat  hydrographic  data-based  estimate  of  Eabrador 
Sea  heat  loss  by  Pickart  and  Spall  (2007)  resulted  in  values  an  order  of  magnitude  smaller. 
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with  the  heat  loss  predominantly  resulting  from  the  horizontal  rather  than  the  overturning 
circulation  both  in  depth  and  density  space.  They  concluded  that  the  Atlantic  MOC  is  not 
largely  impacted  locally  by  deep  convection  in  the  Labrador  Sea.  The  discrepancy  between 
these  results  partly  stems  from  different  locations  and  NAO  periods  from  which  the  hy¬ 
drographic  data  were  obtained,  as  well  as  the  low  horizontal  resolution  of  the  data  used 
by  Talley  (2003).  The  estimate  of  Pickart  and  Spall  (2007)  includes  only  a  limited  part  of 
the  convective  area  in  the  subpolar  gyre  (although  this  does  contain  the  area  where  convec¬ 
tive  activity  is  expected  to  be  most  intense).  In  particular,  convection  in  the  Irminger  Sea 
was  not  accounted  for,  which  may  be  a  considerable  source  of  heat  loss  to  the  atmosphere. 
The  disparity  between  different  observational  estimates  regarding  the  importance  of  LSW 
formation  to  meridional  heat  transport,  not  to  mention  the  disparity  between  observations 
and  GCMs,  illustrates  how  much  remains  to  be  learned  about  the  subpolar  North  Atlantic. 

To  summarize:  The  subpolar  North  Atlantic  is  a  critical  component  of  the  global  cli¬ 
mate  system.  Deep  convection  in  the  subpolar  gyre  releases  heat  from  the  ocean  to  the 
atmosphere,  and  is  crucial  to  the  climate  of  northwestern  Europe  (e.g.  Vellinga  and  Wood, 
2002).  A  resulting  effect  is  that  more  anthropogenic  carbon  is  stored  in  the  subpolar  North 
Atlantic  per  unit  area  than  in  any  other  ocean  (Sabine  et  al,  2004).  Atmospheric  variabil¬ 
ity,  primarily  NAO-related,  strongly  impacts  convective  activity  as  well  as  the  circulation 
and  extent  of  the  subpolar  gyre.  Consequent  changes  in  hydrographic  properties  have  di¬ 
rect  economic  repercussions  on  commercial  fisheries,  even  to  the  extent  that  timeseries  of 
registered  catches  of  various  species  may  be  used  as  proxies  for  the  state  of  the  subpolar 
gyre  (Hatun  et  al,  2009).  Numerical  models  struggle  to  realistically  simulate  the  subpolar 
North  Atlantic  in  various  ways,  but  most  agree  that  LSW  formation  impacts  the  MOC.  The 
Irminger  Sea  is  of  particular  interest  given  its  potential  to  form  LSW. 

Despite  the  evidence  to  date,  there  is  not  an  overall  consensus  in  the  oceanographic 
community  that  formation  of  LSW  can  occur  in  the  Irminger  Sea.  A  second  location  of 
LSW  formation  outside  the  Labrador  Sea  would  have  important  consequences  for  our  un¬ 
derstanding  of  the  North  Atlantic  MOC,  the  modification  of  the  dense  overflow  waters  from 
the  Nordic  seas,  and  the  stratification  and  ventilation  of  the  interior  North  Atlantic.  LSW 
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formation  rates  and  ventilation  times  would  have  to  be  reconsidered,  and  both  observations 
and  model  results  from  this  area  must  be  interpreted  in  a  new  perspective.  Because  of  close 
contact  and  mixing  between  the  NADW  constituent  water  masses,  local  deep  convection, 
and  recipience  of  export  from  the  Arctic,  there  is,  to  quote  Yashayaev  and  Dickson  (2008), 
a  unique  global  importance  of  regional  processes  in  the  Irminger  Sea  for  the  transfer  of 
ocean  climate  signals  between  water  masses  and  to  great  ocean  depths. 

1.6  Thesis  outline 

The  foeus  of  this  thesis  is  on  deep  convection  and  LSW  formation  in  the  subpolar  North 
Atlantic,  with  a  special  emphasis  on  the  Irminger  Sea  and  the  small-scale  atmospheric 
and  oceanic  processes  that  are  important  to  local  overturning.  The  thesis  is  structured  as 
three  self-contained  articles,  two  of  which  are  published;  the  third  will  be  submitted  for 
publication  shortly. 

In  Chapter  2  the  impact  of  Greenland  on  the  atmospheric  circulation  in  its  vicinity  is 
investigated,  with  a  particular  emphasis  on  the  westerly  Greenland  tip  jet.  An  objective  and 
sensitive  method  is  used  to  explore  the  ERA-40  meteorological  reanalysis  product  for  the 
signatures  of  tip  jets,  with  the  goal  of  forming  a  comprehensive  and  consistent  climatology 
of  such  events.  This  is  then  used  to  investigate  the  evolution  of  the  tip  jet  and  to  gain  further 
dynamical  insight  regarding  its  mechanism  of  generation.  The  main  caveat  of  using  a  the 
ERA-40  product  is  its  coarse  resolution,  but  comparison  with  satellite  wind  measurements 
for  the  period  where  both  products  were  available  shows  that  while  the  tip  jet  magnitudes 
were  underestimated,  the  vast  majority  of  the  tip  jets  were  detected.  This  chapter  was  pub¬ 
lished  in  the  Quarterly  Journal  of  the  Royal  Meteorological  Society  as  part  of  the  Greenland 
Elow  Dynamics  Experiment  special  issue.  The  co-authors,  Thomas  Spengler,  Huw  Davies, 
and  Bob  Pickart,  contributed  to  the  study  with  comments  and  insights. 

In  Chapter  3  the  return  of  deep  convection  to  the  subpolar  North  Atlantic  in  winter 
2007-08  following  years  of  shallow  overturning  is  documented.  Various  data  sets  are  ana¬ 
lyzed  with  the  objective  to  elucidate  the  reasons  why  it  happened,  and  the  study  illuminates 
the  complexity  of  the  North  Atlantic  convective  system.  This  chapter  was  published  in  Na- 
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ture  Geoscience.  Specific  contributions  of  the  co-authors  are  as  follows:  Virginie  Thierry 
and  Gilles  Reverdin  first  noticed  the  return  of  deep  convection  from  real-time  ARGO  data 
which  helped  motivate  the  study.  Craig  Lee  and  Brian  Petrie  did  the  freshwater  flux  calcu¬ 
lation  for  Davis  Strait,  which  was  discussed  in  Section  3.6.  Tom  Agnew  and  Amy  Wong 
supplied  movies  of  ice  concentration  and  motion  in  the  Arctic/North  Atlantic  that  moti¬ 
vated  a  part  of  the  discussion  section.  Mads  Ribergaard  provided  meteorological  data  from 
Greenland.  All  of  the  co-authors  contributed  to  the  study  with  comments  and  insights. 

The  main  objectives  of  Chapter  4  are  to  quantify  the  structure,  transport  and  interannual 
variability  of  the  Irminger  Gyre  and  to  investigate  the  importance  of  the  Irminger  Sea  as  a 
source  region  for  LSW.  In  situ  data  were  generously  provided  by  0yvind  Knutsen,  Artem 
Sarafanov,  Herle  Mercier,  Pascale  Lherminier,  Manfred  Bersch,  Hendrik  van  Aken,  Jens 
Meincke,  and  Detlef  Quadfasel. 

Collectively  the  thesis  chapters  are  intended  to  add  to  our  body  of  knowledge  regarding 
the  atmospheric  forcing,  the  nature  and  interplay  of  the  atmosphere-ice-ocean  system,  and 
the  structure  and  variability  of  the  circulation  in  the  subpolar  North  Atlantic,  and  how  these 
phenomena  impact  upon  convection  in  this  region.  The  purpose  is  to  contribute  to  a  basis 
from  which  we  can  better  interpret  existing  observations  and  improve  models,  so  that  we 
can  ultimately  understand  more  thoroughly  how  Earth’s  climate  system  works. 
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Chapter  2 


Multi-event  analysis  of  the  westerly 
Greenland  tip  jet  based  upon  45  winters 
in  ERA-40 


This  chapter  has  been  published  online  as  an  article  in  the  Quarterly  Journal  of  the  Royal 
Meteorological  Society  as  part  of  the  Greenland  Flow  Dynamics  Experiment  special  issue, 
and  is  reprinted  with  their  permission.  The  authors  are  Kjetil  Vage,  Thomas  Spengler,  Huw 
C.  Davies,  and  Robert  S.  Pickart. 

2.1  Abstract 

The  westerly  Greenland  tip  jet  is  an  intense,  narrow  and  intermittent  wind  phenomenon 
located  southeast  of  Cape  Farewell  that  occurs  frequently  during  the  winter  season.  Using 
the  ERA-40  reanalysis  data  set,  a  catalogue  of  586  objectively  detected  westerly  tip  jet 
events  is  compiled  for  the  winters  1957-2002,  and  an  analysis  is  undertaken  of  the  character 
of  the  jet  and  its  accompanying  atmospheric  features.  It  is  shown  that  the  tip  jet  frequency 
exhibits  a  significant  positive  correlation  with  both  the  NAO  index  and  the  latitude  of  the 
Icelandic  Low.  The  peak  wind  speed  and  accompanying  heat  fluxes  of  the  jet  have  values 
up  to  30  m  and  600  W  mT^,  respectively,  and  are  sustained  for  less  than  one  day. 
The  air  parcels  constituting  the  tip  jet  are  shown,  based  upon  a  trajectory  model  and  the 
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ERA-40  data  set,  to  have  a  continental  origin,  and  to  exhibit  a  characteristic  deflection  and 
acceleration  around  southern  Greenland.  The  events  are  almost  invariably  accompanied 
both  by  a  notable  coherence  of  the  lower-level  tip  jet  with  an  overlying  upper-level  jet 
stream,  and  by  a  surface  cyclone  located  to  the  lee  of  Greenland.  It  is  also  shown  that  the 
cyclone  originates  upstream  of  and  is  advected  to  the  lee  of  Greenland,  and  thereby  it  both 
precedes  in  time  and  contributes  dynamically  to  the  formation  of  the  tip  jet.  On  this  basis 
it  is  suggested  that  the  tip  jet  arises  from  the  interplay  of  the  synoptic-scale  flow  evolution 
and  the  perturbing  effects  of  Greenland’s  topography  upon  the  flow. 

2.2  Introduction 

Greenland  exerts  a  significant  impact  on  the  atmospheric  circulation  of  the  Northern  Hemi¬ 
sphere.  Among  other  things,  it  is  known  to  cause  a  southward  shift  of  the  North  Atlantic 
storm  track  (Petersen  et  al,  2004).  Both  the  orography  and  sharp  low-level  temperature 
gradients  arising  from  the  warm  interior  ocean  adjacent  to  the  cold  water  and  ice  on  the 
east  Greenland  shelf  partially  account  for  the  presence  of  the  Icelandic  Low  (Tsukernik 
et  al,  2007).  On  shorter  time  scales,  numerical  simulations  have  shown  that  the  high  to¬ 
pography  of  Greenland  affects  synoptic  storm  systems  in  its  vicinity  (e.g.  Kristjansson  and 
Mclnnes,  1999;  Petersen  et  al,  2003). 

The  region  near  Cape  Farewell  (the  southern  tip  of  Greenland)  is  the  windiest  area  of  the 
World  Ocean  (Sampe  and  Xie,  2007).  This  is  to  a  large  extent  due  to  a  variety  of  intense, 
small-scale  wind  phenomena  that  arise  from  the  impact  of  the  Greenland  topography  on 
low  pressure  systems  that  pass  nearby.  (Such  orographically  induced  winds  are  not  unique 
to  Greenland,  see  for  example  Schwerdtfeger,  1975;  Parish,  1982;  McCauley  and  Sturman, 
1999;  Seefeldt  and  Cassano,  2008).  Near  Greenland  these  phenomena  predominantly  occur 
in  winter.  Most  notable  are  westerly  (forward)  and  easterly  (reverse)  tip  jets  and  barrier 
winds.  Tip  jets  are  narrow,  intermittent  wind  patterns  near  the  southern  tip  of  Greenland, 
while  barrier  winds  are  a  geostrophic  response  to  stable  air  being  forced  towards  the  high 
eastern  coast  of  Greenland  (Moore,  2003;  Moore  and  Renfrew,  2005;  Petersen  et  al,  2009). 
The  focus  of  this  study  is  on  the  westerly  tip  jet  (see  Figure  2- 1 ,  hereafter  referred  to  simply 
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as  a  “tip  jet”). 

Tip  jets  often  develop  when  a  low  pressure  system  moves  into  the  region  east  of  south¬ 
ern  Greenland  (Pickart  et  al,  2003a;  Moore,  2003).  Two  mechanisms  generating  the  tip 
jet  have  been  proposed:  The  first  involves  orographic  descent.  Doyle  and  Shapiro  (1999) 
argued  that  the  tip  jets  are  governed  by  conservation  of  the  Bernoulli  function  when  air 
parcels  descend  sharply  down  the  lee  slope  of  Greenland  (the  Bernoulli  function  is  the  sum 
of  enthalpy,  kinetic  energy,  and  potential  energy  per  unit  mass;  Schar,  1993).  The  second 
mechanism,  also  noted  by  Doyle  and  Shapiro  (1999)  and  later  expanded  upon  by  Moore  and 
Renfrew  (2005),  involves  blocking  due  to  the  high  topography  of  Greenland.  The  dynam¬ 
ical  parameter  space  governing  the  occurrence  of  blocking  in  the  case  when  stratified  flow 
encounters  a  topographic  barrier  has  been  investigated  both  in  the  absence  (Smith,  1989; 
Olafsson  and  Bougeault,  1996)  and  presence  (Pierrehumbert  and  Wyman,  1985;  Petersen 
et  al,  2003)  of  rotation  (see  Doyle  and  Shapiro,  1999,  for  a  theoretical  overview).  Moore 
and  Renfrew  (2005)  found  that  acceleration  through  deflection  of  surface  winds  around  the 
southern  tip  of  Greenland  during  periods  of  blocking  was  a  consistent  feature  of  tip  jets  in 
their  climatological  study  of  high  wind  speed  events  in  the  vicinity  of  Greenland. 

The  small  meridional  scale  of  the  tip  jet  results  in  a  region  of  localized  high  cyclonic 
wind-stress  curl  east  of  Cape  Farewell  (Pickart  et  al,  2003a).  Using  an  idealized  ocean 
model.  Spall  and  Pickart  (2003)  showed  that  this  patch  of  enhanced  wind-stress  curl  is  suffi¬ 
cient  to  drive  the  recirculation  gyre  observed  in  the  Irminger  Sea  east  of  Greenland  (Laven¬ 
der  et  al,  2000).  The  cyclonic  wind-stress  curl  leads  to  upwelling  and  isopycnal  doming 
in  the  water  column  below,  which,  together  with  the  large  heat  fluxes  caused  by  the  strong 
winds,  are  conducive  for  driving  deep  oceanic  convection  in  winter  (Pickart  et  al,  2003a; 
Vage  et  al,  2008).  It  has  been  hypothesized  that  the  southwestern  Irminger  Sea,  within  the 
Irminger  Gyre  and  under  the  direct  influence  of  the  tip  jet,  may  be  another  region  where 
the  densest  type  of  subpolar  mode  water,  Labrador  Sea  Water  (LSW),  is  formed  during  suf¬ 
ficiently  severe  winters  (Pickart  et  al,  2003b).  (LSW  was  previously  thought  to  originate 
solely  from  the  Labrador  Sea).  The  strong  tip  jet  winds  also  present  a  significant  marine 
hazard  (e.g.  Moore,  2003). 

In  this  work  we  use  an  objective  and  sensitive  method  to  explore  the  ERA-40  meteoro- 
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Figure  2-1:  Mean  tip  jet  QuikSCAT  winds  (colour,  m  and  ERA-40  sea  level  pres¬ 
sure  (contours,  hPa)  over  the  period  1999-2002.  See  text  for  explanation  of  rectangle  and 
lines.  The  red  dot  marks  the  location  of  the  Prins  Christian  Sund  (PCS)  meteorological 
station  near  Cape  Farewell.  The  Labrador  Sea  to  the  southwest  and  the  Irminger  Sea  to  the 
southeast  of  Greenland  are  regions  of  dense  water  formation. 


logical  reanalysis  product  for  the  signatures  of  tip  jets.  The  goal  is  to  form  a  comprehensive 
and  consistent  climatology  of  such  events,  which  is  then  used  to  investigate  the  evolution 
of  the  tip  jet  and  to  gain  further  dynamical  insight  regarding  its  mechanism  of  generation. 
The  main  advantages  of  using  a  meteorological  reanalysis  for  climatological  studies  are  the 
long  time  periods  covered  and  the  availability  of  all  relevant  meteorological  fields  on  a  3- 
dimensional  grid,  while  the  main  disadvantage  is  a  relatively  coarse  horizontal  and  vertical 
resolution.  The  advent  of  satellite  scatterometers  (in  particular  the  QuikSCAT  satellite  in 
1999),  which  measure  the  surface  wind  field  over  the  ocean  at  high  resolution,  provides  an 
opportunity  to  check  the  consistency  of  the  lower-resolution  reanalysis  surface  winds  in  the 
vicinity  of  Cape  Farewell. 

The  paper  is  outlined  as  follows.  The  ERA-40  reanalysis  product  and  the  tip  jet  de¬ 
tection  routine  are  described  in  Section  2.3.  A  comparison  with  the  QuikSCAT  satellite 
scatterometer  winds  follows  in  Section  2.4.  In  Section  2.5  the  relationship  between  tip  jet 
frequency  and  large-scale  North  Atlantic  atmospheric  patterns  are  examined.  Composite 


23 


averages  provide  a  robust  description  of  the  evolution  of  the  tip  jet  and  are  described  in 
Section  2.6.  A  Lagrangian  trajectory  analysis,  presented  in  Section  2.7,  complements  the 
description  of  the  tip  jet  by  identifying  the  paths  and  origin  of  the  air  parcels  contained 
within  the  jet.  The  atmospheric  conditions  during  tip  jet  events  are  investigated  in  Sec¬ 
tion  2.8.  In  Section  2.9  we  discuss  the  dependency  of  the  tip  jet  on  the  low-level  flow 
interaction  with  the  topography  of  Greenland  as  well  as  on  the  large-scale  atmospheric 
conditions. 


2.3  Data  and  methods 

2.3.1  ERA-40  reanalysis  product 

The  European  Centre  for  Medium-Range  Weather  Forecasts  (ECMWF)  45-year  reanalysis 
product  (ERA-40,  Uppala  et  al,  2005)  covering  the  period  1957  to  2002  was  analyzed  in 
this  study.  The  horizontal  and  temporal  resolutions  at  the  surface  and  at  23  pressure  lev¬ 
els  are  one  degree  and  six  hours,  respectively.  Produced  with  a  fixed  numerical  weather 
prediction  data  assimilation  system,  ERA-40  (and  other  reanalysis  products)  provides  a 
physically  consistent  interpolation  of  historical  observational  data  through  space  and  time. 
However,  the  observing  system  did  not  remain  constant  during  this  period.  The  growing 
importance  of  satellite  instruments  after  1979  had  a  particularly  large  impact  on  the  meteo¬ 
rological  analyses,  mainly  in  regions  where  ground-based  observations  were  sparse  and  in 
the  upper  troposphere  and  stratosphere  (Bengtsson  et  al,  2004).  Thus  the  ERA-40  fields 
may  not  be  regarded  as  fully  consistent  throughout  the  entire  period.  For  our  tip  jet  cli¬ 
matology  this  issue  is  not  addressed,  but  it  is  less  of  a  concern  in  the  northern  hemisphere. 
Furthermore  upper-level  atmospheric  data  are  only  considered  for  the  most  recent  period. 

As  a  check  on  the  accuracy  of  the  ECMWF  model  in  a  subpolar  region  of  sparse 
data,  Renfrew  et  al.  (2002)  compared  the  surface  fields  from  the  National  Centers  for  En¬ 
vironmental  Prediction  (NCEP,  Kalnay  et  al,  1996)  and  the  ECMWF  operational  analysis 
with  data  recorded  by  the  research  vessel  Knorr  during  a  1997  winter  cruise  in  the  Eabrador 
Sea.  They  found  that  the  ECMWF  model  fields  represented  the  observed  data  reasonably 
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well,  but  contained  a  cold  bias  for  sea  surface  temperature  and  2  m  air  temperature  near 
the  ice  edge,  and  a  tendency  to  overestimate  the  latent  and  sensible  heat  fluxes.  In  the 
present  study  the  latter  issue  is  circumvented  to  some  extent  by  application  of  a  bulk  for¬ 
mula  (Fairall  et  ai,  2003)  to  compute  the  heat  fluxes  in  Section  2.6. 

Concerns  regarding  the  ability  of  a  global  product  of  relatively  coarse  resolution  like 
ERA-40  to  accurately  capture  small-scale  atmospheric  phenomena  such  as  the  Greenland 
tip  jet  (illustrated  in  Figure  2-2)  have  been  raised  previously  (Pickart  et  al.,  2003a;  Moore, 
2003).  This  is  examined  further  in  Section  2.4. 

QuikSCAT  winds  ERA-40  winds 
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Figure  2-2:  Surface  wind  field  during  the  February  21,  2002,  tip  jet  event  illustrat¬ 

ing  the  difference  between  (a)  QuikSCAT  and  (b)  ERA-40. 

2.3.2  Tip  jet  detection 

Two  empirical  orthogonal  function  (EOF)  approaches  have  been  used  in  previous  studies  to 
objectively  detect  tip  jet  events.  Pickart  et  al.  (2003a)  found  that  tip  jets  were  characterized 
by  strong  westerly  winds,  anomalously  low  air  temperatures,  and  low  sea  level  pressure. 
They  constructed  an  EOF  based  on  these  three  variables  using  nearly  30  years  of  meteo¬ 
rological  data  from  the  Prins  Christian  Sund  weather  station  near  Cape  Farewell  (labeled 
PCS  in  Fig.  2-1)  to  identify  the  events.  Vage  et  al.  (2008)  noted  that  tip  jet  events  were 
relatively  insensitive  to  the  location  and  central  pressure  of  their  parent  low  pressure  sys¬ 
tems,  and  showed  that  the  gradient  of  sea  level  pressure  was  a  better  metric  than  sea  level 
pressure  alone.  Some  tip  jets  were  also  observed  to  completely  evade  the  meteorological 
station.  Correspondingly,  Vage  et  al.  (2008)  devised  an  EOF  routine  that  employed  time- 
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series  of  (i)  the  maximum  10  m  zonal  wind  speed  within  the  blue  rectangle  of  Figure  2-1 
recorded  by  the  QuikSCAT  satellite  scatterometer  (retrieved  using  the  “Ku-2001”  model, 
Wentz  et  al,  2001);  (ii)  the  mean  NCEP  sea  level  pressure  gradient  along  the  lines  a,  b, 
and  c  in  Figure  2-1;  and  (iii)  the  PCS  air  temperature.  This  method  proved  to  be  more 
sensitive  to  tip  jet  detection,  and  had  the  advantage  of  utilizing  accurate,  high-resolution 
satellite  wind  speed  data  covering  a  larger  geographical  area.  Vage  et  al.  (2008)  used  a 
threshold  zonal  wind  speed  of  25  m  5  ^  to  define  a  robust  tip  jet.  Manual  inspection  of  the 
scatterometer  winds  confirmed  that  all  of  the  robust  tip  jets  that  took  place  during  the  two 
winters  that  they  considered  were  detected. 

In  this  work  we  primarily  use  the  EOF  routine  developed  by  Vage  et  al.  (2008).  We  note 
that  satellite  scatterometers  were  not  available  prior  to  1987,  and  that  there  are  temporal 
gaps  in  the  PCS  air  temperature  record.  Hence  we  use  only  timeseries  extracted  from 
ERA-40  as  input  for  the  EOE  routine.  ERA-40  was  used  alone  in  order  to  have  a  consistent 
detection  algorithm  covering  the  entire  period.  Eurthermore  it  is  desirable  to  avoid  data 
from  the  location  of  the  PCS  weather  station  due  to  suspected  influence  of  the  high  local 
topography.  Instead  of  temperatures  from  the  PCS  location,  we  used  a  roving  temperature 
timeseries.  At  each  time  step,  the  temperature  at  the  location  of  the  maximum  zonal  wind 
speed  within  the  blue  rectangle  of  Eigure  2- 1  was  tabulated.  In  this  way  the  temperature 
at  the  core  of  the  tip  jet  was  always  considered.  Eor  completeness  we  also  applied  the 
same  EOE  method  implemented  by  Pickart  et  al.  (2003a)  using  the  PCS  data  only.  Eor  the 
majority  of  cases  the  two  EOE  routines  identified  the  same  events. 

While  we  use  the  same  routine  as  Vage  et  al.  (2008),  our  reconstructed  tip  jet  winds  are 
weaker  than  theirs  largely  because  of  the  coarser  resolution  of  the  ERA-40  winds.  As  such, 
the  threshold  reconstructed  zonal  wind  speed  for  defining  an  event  as  robust  was  lowered 
from  25  to  18  m  .  This  limit  was  chosen  to  achieve  a  balance  between  avoiding  false 
positive  identifications,  and  at  the  same  time  not  missing  robust  tip  jets.  The  threshold 
value  was  determined  by  carefully  comparing  the  ERA-40  and  QuikSCAT  winds  for  the 
period  of  overlap  (1999-2002).  In  reality  there  exists  a  continuum  of  tip  jet  magnitudes, 
and  the  label  “robust”  is  somewhat  arbitrary.  All  tip  jets  discussed  from  this  point  onwards 
are  robust  according  to  the  above  definition,  hence  the  label  will  be  dropped. 
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All  of  the  events  that  satisfied  the  criteria  from  either  of  the  two  EOF  routines  above 
were  examined  manually  to  ensure  that  each  tip  jet  was  recorded  once,  at  its  peak  intensity, 
and  to  locate  the  centres  of  their  parent  low  pressure  systems.  In  addition  to  the  two  EOFs,  a 
final  procedure  was  applied  to  detect  wind  speed  events  exceeding  20  with  a  direction 
deviating  at  most  30°  from  west  within  the  blue  rectangle  of  Figure  2-1.  This  was  done  to 
identify  any  events  that  might  have  been  missed  by  the  EOFs.  Only  a  few  such  cases  were 
identified.  In  total,  586  tip  jet  events  were  detected  in  the  ERA-40  reanalysis  for  the  winters 
(Nov-Apr)  of  1957  to  2002.  This  is  on  average  13  per  winter.  Validation  of  our  climatology 
through  comparison  with  QuikSCAT  for  the  period  of  overlap  indicates  that  88%  of  all  of 
the  tip  jets  that  had  mean  QuikSCAT  wind  speeds  within  the  blue  rectangle  in  Figure  2-1 
exceeding  25  m  were  detected  in  ERA-40  without  reporting  false  positives.  The  tip 
jets  found  only  in  the  QuikSCAT  data  set  were  in  general  weaker  and/or  of  shorter  duration 
than  the  tip  jets  identified  in  ERA-40. 

2.4  ERA-40  versus  QuikSCAT  winds 

2.4.1  Winter  conditions 

To  examine  the  correspondence  between  ERA-40  and  QuikSCAT  winds,  the  mean  10  m 
wind  speed  within  the  blue  rectangle  of  Figure  2- 1  was  determined  for  every  realization 
of  each  data  set  during  the  winters  of  overlap  (1999-2002).  The  ERA-40  timeseries  was 
then  temporally  interpolated  onto  the  same  twice  daily  time  base  as  the  QuikSCAT  data 
(Fig.  2-3a).  The  correlation  between  the  wind  speeds  from  the  two  products  in  this  region 
is  0.89.  All  correlations  reported  are  significant  at  the  99%  confidence  level.  (Confidence 
intervals  were  determined  using  a  bootstrap  algorithm,  a  procedure  that  involves  random 
sampling  with  replacement  from  the  data  set  and  does  not  require  any  assumptions  about 
the  underlying  probability  distribution.)  The  linear  least  squares  best  fit  resulted  in  the 
following  relation  between  the  two  wind  speed  products: 

\^\QuikSCAT  =  l-4|v|£iM-40-0-9  (2.1) 
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Figure  2-3:  Relationships  between  (a)  mean  ERA-40  and  QuikSCAT  10  m  wind  speeds 
within  the  blue  rectangle  of  Figure  2- 1  twice  daily  during  all  of  the  winters  of  overlap  and 
(b)  maximum  ERA-40  and  QuikSCAT  10  m  wind  speeds  found  within  the  same  box  during 
tip  jet  events  (43  events  in  total). 


In  the  region  near  Cape  Farewell  the  ERA-40  and  QuikSCAT  winds  are  generally  well 
correlated,  but  the  good  correlation  tends  to  break  down  for  the  highest  wind  speeds.  In 
particular,  all  of  the  QuikSCAT  wind  speeds  exceeding  30  m  were  associated  with 
small-scale  wind  events:  Westerly  and  easterly  tip  jets  and  barrier  winds.  As  mentioned 
in  the  introduction,  these  wind  phenomena  result  from  the  impact  of  the  high  topography 
of  Greenland  on  a  nearby  low  pressure  system,  and  it  is  clear  that  they  are  not  as  well 
represented  in  the  lower  resolution  ERA-40  fields.  However,  Moore  et  al.  (2008)  found 
that,  relative  to  a  meteorological  buoy  in  the  Irminger  Sea,  QuikSCAT  showed  evidence  of 
a  high  wind  speed  bias.  Renfrew  et  al.  (2009)  reached  a  similar  conclusion  using  aircraft- 
based  observations  east  of  southern  Greenland.  This  bias  may  be  a  contributing  factor  to 
the  reduced  correlation  at  extreme  wind  speeds,  though  the  underestimation  of  high  winds 
in  ERA-40  due  to  low  resolution  and  orographic  effects  is  likely  the  dominant  reason. 
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2.4.2  Tip  jet  conditions 

Considering  only  the  wind  speeds  during  tip  jet  events,  the  correlation  between  ERA-40 
and  QuikSCAT  maximum  wind  speeds  dropped  to  0.75  (Fig.  2-3b),  and  the  least  squares 
relation  changed  to 


\y\QuikSCAT  =  3.1|v|£iM-40  —  30  (2.2) 

We  chose  to  consider  the  maximum  wind  speed  instead  of  the  mean  because  the  scale  of 
the  tip  jet  is  much  smaller  than  the  rectangle.  (When  considering  the  mean  wind  speed 
during  tip  jet  events,  the  correlation  was  0.78.)  We  will  also  use  (2.2)  to  estimate  the  “true” 
wind  speed  at  the  core  of  the  composite  average  tip  jet  in  Section  2.6. 

It  is  clear  that  ERA-40  significantly  underestimates  the  wind  speeds  during  tip  jet  events 
relative  to  the  scatterometer.  Using  the  global  ECMWF  forecast  model,  Jung  and  Rhines 
(2007)  found  that  increasing  the  model  resolution  beyond  that  of  ERA-40  continuously 
improved  the  representation  of  the  tip  jet  as  judged  by  high-resolution  scatterometer  studies 
of  the  tip  jet  (Moore  and  Renfrew,  2005).  In  particular,  both  the  heat  flux  and  surface  wind 
stress  curl  were  increased.  In  any  case,  the  significant  positive  correlation  between  ERA-40 
and  QuikSCAT  winds  during  tip  jet  events  indicates  that  our  tip  jet  climatology  effectively 
captures  the  most  intense  tip  jet  events. 

2.5  Tip  jets  and  the  NAO 

The  North  Atlantic  Oscillation  (NAO)  is  the  dominant  pattern  of  climate  variability  in  the 
North  Atlantic  (Hurrell,  1995).  Winters  that  are  characterized  by  a  positive  NAO  index  are 
associated  with  stronger  westerly  winds,  enhanced  air-sea  buoyancy  fluxes,  and  are  often 
linked  with  deep  convection  in  the  subpolar  North  Atlantic  Ocean  (Dickson  et  al,  1996; 
Vage  et  al,  2009a).  Pickart  et  al  (2003a)  also  found  that  a  greater  number  of  tip  jet  events 
took  place  during  high  NAO  index  winters.  Figure  2-4  shows  the  timeseries  of  winter 
(Nov-Apr)  tip  jet  frequency  over  the  45  year  period  of  ERA-40.  Significant  interannual 
variability  is  evident;  some  winters  have  as  few  as  3-5  tip  jets,  while  others  have  as  many 
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as  20-25.  On  average,  there  are  13  tip  jet  events  per  winter,  with  a  standard  deviation  of  5. 
Figure  2-4  also  illustrates  the  temporal  relation  between  tip  jet  frequency  and  the  monthly 
NAO  index  of  Hurrell  (1995)  averaged  over  the  winter  months  (Nov- Apr).  Using  the  tip  jet 
climatology  of  Pickart  et  al.  (2003a)  and  the  atmospheric  NAO  pattern  decomposed  into 
indices  of  Icelandic  Low  (IL)  and  Azores  High  (AH)  centres  of  actions  (COA,  Hameed 
et  al,  1995),  Bakalian  et  al.  (2007)  found  that  the  number  of  tip  jet  events  per  winter  was 
more  sensitive  to  the  latitude  of  the  IL  than  to  the  NAO  index. 


Using  the  present  climatology,  with  the  monthly  NAO  index  of  Hurrell  (1995)  and  the 
COA  indices  of  Hameed  et  al.  (1995)  averaged  over  the  winter  months  (Nov-Apr),  we  find 
similar  correlations.  This  is  true  despite  the  underrepresented  tip  jet  frequency  in  Pickart 
et  al.  (2003a) ’s  climatology  (see  Vage  et  al.,  2008).  In  particular,  we  find  correlations  of 
0.71  and  0.69  between  the  number  of  tip  jet  events  per  winter  and  the  NAO  index  and  lat¬ 
itude  of  the  IL,  respectively.  This  suggests  that  the  NAO  state  and  the  meridional  location 
of  the  IL  are  of  comparable  importance  for  the  frequency  of  tip  jet  events  (Fig.  2-5).  As  the 
location  of  the  IL  may  influence  the  value  of  the  station-based  NAO  index,  these  variables 
are  not  expected  to  be  completely  independent.  Indeed,  a  correlation  of  0.86  was  found  be- 
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tween  the  NAO  index  and  the  latitude  of  the  IL.  Using  the  wintertime  NAO  index  (Hurrell, 
1995)  based  upon  the  months  Deeember  through  Mareh,  the  eorrelation  decreased  to  0.64 
(a  similarly  low  correlation  of  0.65  was  obtained  when  averaging  the  monthly  NAO  index 
over  the  same  months).  This  shows  that  the  shoulder  months,  November  and  April,  which 
are  important  from  an  oceanic  convection  perspective,  also  contribute  a  sizeable  number  of 
tip  jet  events  (16  and  8  percent,  respectively). 


NAO  index  IL  latitude 

Figure  2-5:  Relationships  between  winter  tip  jet  frequency  and  (a)  the  NAO  index  and  (b) 
the  latitude  of  the  Icelandic  Low. 


2.6  Tip  jet  composites 

Composite  averages  of  all  586  tip  jets  from  the  ERA-40  data  set  during  the  24  hour  period 
surrounding  an  event  portray  the  tip  jet  as  an  intense,  short-lived  phenomenon  (Fig.  2- 
6).  Peak  wind  speeds  approaching  20  m  were  sustained  for  less  than  a  day  (Fig.  2- 
6a).  Pickart  et  al.  (2003a)  found  that  tip  jets  at  the  Prins  Christian  Sund  meteorological 
station  on  average  lasted  approximately  3  days.  The  discrepancy  between  these  estimates 
of  tip  jet  duration  is  mainly  due  to  our  focus  on  peak  tip  jet  wind  speeds,  while  they  also 
considered  the  shoulders  of  the  tip  jet  events.  Using  Equation  2.2  the  “true”  wind  speed 
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at  the  core  of  the  composite  average  tip  jet  was  estimated  to  exceed  30  m  The  region 
of  low  sea  level  pressure  localized  northeast  of  Cape  Farewell  was  at  its  deepest  and  most 
compact  at  the  peak  of  the  tip  jets  (Fig.  2-6b),  and  at  a  lag  of  +12  hours  had  moved  in  the 
downstream  direction  towards  Iceland.  This  is  reflected  in  the  gradient  of  sea  level  pressure. 
The  strong  pressure  gradient  extending  eastward  from  Cape  Farewell,  which  in  agreement 
with  Vage  et  al.  (2008)  was  found  to  be  a  feature  common  to  every  tip  jet  event,  was 
most  intense  at  the  peak  of  the  composite  event  compared  to  lags  plus  and  minus  12  hours. 
This  is  consistent  with  the  composite  velocities  (Fig.  2-6a).  In  agreement  with  previous 
studies  (e.g.  Pickart  et  al,  2003a),  every  tip  jet  event  was  associated  with  a  parent  low 
pressure  system  (marked  with  an  L  in  Figure  2-6b).  Although  most  of  the  parent  cyclones 
were  located  close  to  the  east  coast  of  Greenland,  there  was  a  significant  scatter  in  their 
positions,  but  never  to  the  extent  that  the  strong  pressure  gradients  east  of  Cape  Farewell 
were  compromised. 

Turbulent  heat  fluxes  were  computed  using  a  bulk  formula  (COARE  3.0,  Fairall  et  al, 
2003)  with  the  inputs  of  wind  speed,  humidity,  air  temperature,  and  sea  surface  temper¬ 
ature,  and  were  on  average  more  than  3  times  greater  during  tip  jet  events  compared  to 
background  levels  (Fig.  2-6c).  Replacing  the  ERA-40  wind  speed  with  our  estimate  of  the 
true  wind  speed  in  the  bulk  formula  raised  the  heat  fluxes  in  the  centre  of  the  composite 
tip  jet  from  400  W  rrT^  to  almost  600  W  mT^.  Even  though  the  QuikSCAT-based  true 
wind  speed  estimate  may  have  a  positive  bias  (Moore  et  al,  2008;  Renfrew  et  al,  2009), 
the  magnitude  of  the  heat  flux  at  the  core  of  the  tip  jet  is  still  likely  to  be  underestimated. 
Tip  jets  are  associated  with  a  drop  in  air  temperature  (Pickart  et  al,  2003a),  and  it  will  be 
shown  in  Section  2.7  that  cold,  low-humidity  air  of  continental  origin  tends  to  comprise  the 
tip  jets.  Neither  of  these  features  are  well  captured  by  the  ERA-40  data  set,  and  both  would 
tend  to  further  increase  the  tip  jet  heat  fluxes.  This  suspected  heat  flux  underestimate  is 
also  consistent  with  the  case  studies  of  Jung  and  Rhines  (2007). 

As  seen  in  Eigure  2-6c,  large  heat  fluxes  are  also  present  over  the  Eabrador  Sea  12 
hours  prior  to  the  peak  of  the  tip  jet  events.  These  are  most  likely  associated  with  cold 
air  outbreaks,  during  which  low  pressure  systems  following  the  North  Atlantic  storm  track 
draw  cold  and  dry  continental  air  from  Eabrador  over  the  relatively  warm  ocean,  thereby 
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Figure  2-6:  Composite  averages  of  (a)  wind  speed  (m  5^^),  (b)  sea  level  pressure  (hPa), 
and  (c)  turbulent  heat  fluxes  (VF  m  illustrating  the  temporal  evolution  of  the  Greenland 
tip  jet.  The  title  of  each  panel  refers  to  the  time  lag  relative  to  the  peak  of  the  composite  tip 
jet.  The  L’s  in  (b)  indicate  the  locations  of  the  centres  of  the  parent  cyclones,  the  sizes  of 
the  L’s  are  proportional  to  the  number  of  cyclones  that  they  represent. 
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enhancing  the  air-sea  heat  fluxes  (e.g.  Pickart  et  al,  2008).  The  analysis  of  Vage  et  al. 
(2009a)  shows  that  the  very  same  storms  that  cause  cold  air  outbreaks  in  the  Labrador  Sea 
tend  to  trigger  tip  jet  events  a  few  hours  later  as  they  progress  northeastward  along  the 
storm  track.  Both  of  these  locations  of  high  heat  fluxes,  the  western  Labrador  Sea  and 
the  southwestern  Irminger  Sea,  are  regions  where  deep  convection  is  known  to  take  place 
during  severe  winters.  For  the  evolution  of  the  oceanic  mixed  layer  in  the  Irminger  Sea,  the 
integrated  effect  of  many  tip  jets  during  the  course  of  a  winter  is  of  great  importance  (Vage 
et  al,  2008). 

2.7  Air  parcel  trajectories 

To  better  understand  the  potential  impact  of  tip  jets  on  oceanic  convection,  it  is  necessary 
to  know  the  origin  of  the  air  that  they  transport.  Cold  and  dry  continental  air  will  extract 
more  heat  from  the  ocean  than  relatively  warm  and  moist  maritime  air.  The  upstream  path 
of  air  parcels  comprising  the  tip  jets  might  also  provide  evidence  supporting  one  of  the 
proposed  mechanisms  of  formation;  i.e.  either  the  orographic  descent  hypothesis  of  Doyle 
and  Shapiro  (1999)  or  the  acceleration  through  deflection  hypothesis  advocated  by  Moore 
and  Renfrew  (2005). 

A  3-dimensional  Lagrangian  trajectory  model  (Wernli  and  Davies,  1997)  was  applied 
to  compute  2-day  backwards  air  parcel  trajectories  terminating  above  the  southwestern 
Irminger  Sea  (900  to  950  hPa)  using  ERA-40  winds  as  input.  Only  trajectories  from  the 
core  of  each  tip  jet  were  included,  selected  objectively  by  requiring  the  wind  speed  to 
exceed  20  m  5  ^  and  the  direction  to  deviate  no  more  than  30°  from  west  at  the  point  of 
termination.  Given  the  evolving  observing  system  (Bengtsson  et  al,  2004)  and  computation 
constraints,  only  the  most  recent  winters  (1994  to  2002)  were  considered  for  this  analysis, 
resulting  in  nearly  3000  trajectories  from  101  tip  jet  events  (Fig.  2-7). 

The  trajectories  so  computed  clearly  illustrate  the  continental  origin  of  the  air  parcels 
constituting  the  tip  jet.  This  has  important  implications  for  the  air-sea  heat  exchange  in  the 
Irminger  Sea,  as  the  air  originating  from  northern  Canada  in  winter  is  usually  very  cold 
and  dry.  Before  reaching  the  Irminger  Sea,  the  air  must  cross  the  Labrador  Sea,  whose 
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Figure  2-7:  Backward  2-day  trajectories  of  air  parcels  within  tip  jets  (see  text  for  selec¬ 
tion  parameters).  Note  the  continental  origin  of  the  air  composing  the  tip  jet  and  how  the 
trajectories  curve  around  southern  Greenland. 

western  margin  and  northern  end  are  covered  with  ice  in  winter.  Hence  modification  by 
air-sea  interaction  in  the  Labrador  Sea  is  somewhat  reduced,  and,  as  shown  in  Section  2.6, 
tip  jets  are  very  effective  at  extracting  heat  from  the  southwestern  part  of  the  Irminger  Sea. 
A  marked  difference  in  origin  of  the  air  comprising  the  westerly  and  easterly  tip  jets  is 
likely  a  major  reason  why  the  latter  is  not  believed  to  cause  deep  oceanic  convection  in  the 
southeastern  Labrador  Sea  (Sproson  et  al,  2008). 

Figure  2-7  also  demonstrates  that  most  of  the  trajectories  curve  around  Greenland,  sug¬ 
gesting  that  acceleration  by  deflection  is  a  key  mechanism  of  tip  jet  generation  (e.g.  Moore 
and  Renfrew,  2005).  In  order  to  shed  more  light  on  this,  the  pressure  change  and  veloc¬ 
ity  along  each  trajectory  were  computed  and  interpolated  onto  a  0.5°  x  0.5°  grid  using  a 
Laplacian- spline  interpolator.  The  resulting  along-trajectory  pressure  change  (Fig.  2-8a, 
colour)  indicates  that  some  sinking  takes  place  over  southern  Greenland,  which  would  sup¬ 
port  the  orographic  descent  under  Bernoulli  conservation  hypothesis  of  Doyle  and  Shapiro 
(1999).  However,  this  involves  only  a  small  number  of  trajectories  (less  than  18%  of  the 
trajectories  experienced  a  temporal  pressure  change  greater  than  25  hPa  hr~^.  Fig.  2-8, 
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contours).  The  vast  majority  of  the  trajectories  pass  to  the  south  of  Cape  Farewell,  and  the 
acceleration  associated  with  deflection  around  Greenland  is  evident  (Fig.  2-8b,  colour). 
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Figure  2-8:  Pressure  change  (a,  hPa  hr~^)  and  velocity  (b,  m  5'”^)  computed  along  the 
trajectories  from  Figure  2-7.  The  contours  denote  the  number  of  realizations,  and  the  black 
line  marks  the  mean  tip  jet  trajectory. 


The  mean  Froude  number  estimated  west  of  Greenland  during  tip  jet  events  is  much  less 
than  1.  (The  Froude  number  is  computed  as  Fr  =  where  U  is  wind  speed,  N  is  buoy¬ 
ancy  frequency,  a  measure  of  static  stability,  and  H  is  the  obstacle  height.  U  and  N  were 
computed  from  the  ERA-40  climatology  upstream  and  below  the  height  of  the  obstacle, 
which  was  taken  to  be  2500  m.)  This  indicates  that  stratification  dominates  the  inertia  of 
the  air  parcels,  and  hence  the  air  is  more  likely  to  flow  around  than  over  Greenland  (Skeie 
et  ah,  2006).  This  corroborates  the  evidence  from  the  trajectory  model. 
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2.8  Flow  features  accompanying  the  tip  jet 


To  improve  our  understanding  of  the  processes  leading  to  the  formation  of  tip  jets,  we 
examine  the  major  atmospheric  flow  features  that  accompany  these  events  and  how  they 
relate  to  the  results  from  the  air  parcel  trajectory  model.  In  particular  we  focus  on  the 
low-level  cyclones  and  the  tropopause-level  jet  stream. 

2.8.1  Low-level  cyclone 

It  was  noted  earlier  that  the  synoptic-scale  flow  incident  upon  Greenland  can  generate  nu¬ 
merous  small-scale  weather  phenomena  (see  Renfrew  et  al,  2008,  for  an  overview).  Here 
we  examine  the  occurrence  and  nature  of  the  low-level  cyclonic  systems  that  appear  to  be  a 
ubiquitous  accompanying  feature  of  tip  jets.  The  aim  is  to  shed  light  upon  the  relationship 
between  the  cyclone  and  the  tip  jet.  To  this  end  we  use  a  data  base  (Wernli  and  Schwierz, 
2006)  that  identifies  the  storm  tracks  and  cyclone  field  (defined  as  the  area  within  the  out¬ 
ermost  closed  2  hPa  pressure  contour  of  a  low  pressure  system)  based  upon  the  ERA-40 
reanalysis  sea  level  pressure  data.  To  select  the  cyclones  that  were  associated  with  tip  jet 
events  from  the  data  base,  we  matched  their  spatio-temporal  locations  with  the  manually 
determined  positions  of  the  tip  jets’  parent  low  pressure  systems.  A  success  rate  of  98% 
testifies  to  the  robustness  of  the  data  base. 

Consider  in  turn  the  storm  track  and  the  along-track  pressure  change  and  speed  of  trans¬ 
lation  of  the  subset  of  ERA-40  cyclones  that  were  associated  with  tip  jet  events.  Eor  the 
tracks  (Eig.  2-9)  it  is  evident  that  most  of  the  cyclones  linked  to  the  occurrence  of  tip  jets 
followed  a  well-defined  path  from  the  eastern  seaboard  of  North  America  into  the  Irminger 
Sea.  A  similar  picture  emerged  from  the  manually  tracked  parent  low  pressure  systems  for 
the  two  individual  winters  of  2002-03  and  2003-04  analyzed  by  Vage  et  al.  (2008),  and  is 
also  consistent  with  the  results  of  Hoskins  and  Hodges  (2002). 

The  pre-existence  of  the  cyclonic  systems  prior  to  their  arrival  over  the  Irminger  Sea  has 
important  ramifications.  It  demonstrates  that  the  cyclogenesis  that  spawns  these  systems 
tends  to  occur  upstream  of  Greenland  rather  than  being  an  example  of  in-situ  lee  cycloge¬ 
nesis  off  southern  Greenland  (e.g.  Jung  and  Rhines,  2007;  Kristjansson  and  Thorsteinsson, 
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Figure  2-9:  Tracks  of  tip  jet-generating  cyclones  for  the  duration  of  their  existence.  The 
L’s  mark  the  locations  of  the  centres  of  the  parent  cyclones  at  the  peak  of  each  event,  with 
the  size  proportional  to  the  number  of  cyclones  represented. 

2009;  Mclnnes  et  al,  2009).  Consequently  it  suggests  that  orography  alone  does  not  pro¬ 
duce  the  conditions  necessary  for  tip  jet  formation,  and  that  the  upstream  synoptic-scale 
flow  evolution  plays  an  important  role  as  well. 

To  examine  these  issues  further,  consider  the  pressure  change  and  translation  speed  of 
these  systems  (Fig.  2-10)  computed  along  each  cyclone’s  track  in  the  same  manner  as  the 
along-trajectory  evolution  in  Section  2.7.  These  patterns  show  that  the  mean  behaviour  of 
the  cyclones  was  to  deepen  while  rapidly  approaching  Greenland  from  the  southwest,  and 
thereafter  decelerate  and  start  to  fill  as  they  traverse  the  Irminger  Sea  and  migrate  toward 
Iceland.  The  large  translation  speeds  northwest  of  Iceland  in  Figure  2- 10b  reflects  a  ten¬ 
dency  of  a  substantial  subset  of  cyclones  to  rapidly  transit  through  the  Denmark  Strait  sep¬ 
arating  Greenland  and  Iceland.  These  features  both  indicate  that  the  storms  were  spawned 
upstream  and  de-emphasize  the  role  of  lee-cyclogenesis,  either  acting  in  isolation  or  in  tan¬ 
dem  with  the  synoptic  development,  on  the  subset  of  low  pressure  systems  associated  with 
tip  jet  events.  Also,  the  decrease  in  translation  speed  of  the  systems  between  Greenland 
and  Iceland  (Fig.  2- 10b  and  Hoskins  and  Hodges,  2002)  is  consistent  with  this,  suggest- 
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Figure  4-6:  Satellite  altimeter  ground  tracks  provided  by  AVISO  and  plotted  in  Google  Earth.  The 
red  tracks  belong  to  Topex/Poseidon  and  Jason  (9.1  day  repeat  cycle),  the  yellow  tracks  to  ERS  and 
Envisat  (35  days  repeat  cycle),  and  the  green  tracks  to  GEO  (17  days  repeat  cycle). 

2006).  From  the  monthly  mean  winds  we  computed  wind  stress  curl  using  the  method 
of  Milliff  and  Morzel  (2001).  The  spatial  resolution  of  the  NARR  fields  is  32  km.  The 
increased  resolution  relative  to  the  regular  global  reanalysis  products  (e.g.  NCEP)  is  im¬ 
portant  for  a  more  realistic  rendition  of  the  tip  jet  (Jung  and  Rhines,  2007). 

4.3.5  Projection  and  gridding 

In  order  to  facilitate  the  analysis  and  enable  the  use  of  quantitative  techniques  such  as  em¬ 
pirical  orthogonal  functions  (EOFs),  the  hydrographic  sections  were  first  interpolated  onto 
a  standard  grid  in  the  x:  —  z  plane.  This  was  accomplished  as  follows.  A  least  squares  regres¬ 
sion  line  was  fit  through  the  locus  of  the  hydrographic  station  positions  (stations  from  the 
sections  along  59.5°N  and  60°N  were  excepted).  As  the  majority  of  the  sections  followed 
the  WOCE  A1E/AR7E  line  closely,  the  regression  line  is  nearly  identical  to  the  WOCE 
line,  and  the  regression  line  will  hereafter  be  referred  to  simply  as  AR7E.  A  bathymetry 
profile  along  AR7E  was  extracted  from  Etopo2  (Eig.  4-1).  A  second  bathymetry  profile 
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was  estimated  from  the  maximum  depth  of  each  station  under  the  assumption  that  the 
measurements  ended  10  m  above  the  bottom  (e.g.  Pickart,  2004).  Overall  there  was  very 
good  agreement  between  the  two  resulting  depth  curves.  The  final  bathymetry  profile  is 
a  smoothed  combination  of  the  Etopo2  and  station-based  depth  estimates,  with  the  deeper 
profile  prevailing  at  each  distance  interval  (Fig.  4-7). 

Every  station  was  subsequently  projected  onto  the  AR7E  line.  To  the  west  of  the  center 
of  the  Irminger  Basin  the  location  of  each  station  along  AR7E  was  tentatively  determined 
by  minimizing  the  distance  between  its  original  position  and  the  AR7E  line.  In  the  cases 
where  the  discrepancy  between  the  station’s  original^  and  projected  depths  exceeded  a 
set  tolerance  criterion  (5%),  the  projected  location  was  moved  along  the  AR7E  line  until 
its  original  depth  matched  the  bathymetry  profile  within  the  set  tolerance.  In  the  eastern 
part  of  the  Irminger  Basin  the  northward  shoaling  of  the  Reykjanes  Ridge  towards  Iceland 
presented  a  complication.  As  the  topography  of  the  Reykjanes  Ridge  has  an  important  dy¬ 
namical  effect  on  the  flow  (see  for  instance  Fig.  4-5),  it  was  decided  that  the  distance  of 
each  station  from  the  center  of  the  Reykjanes  Ridge  must  be  maintained  under  the  projec¬ 
tion.  Thus  in  the  eastern  part  of  the  Irminger  Basin  each  station  was  first  projected  onto 
the  AR7E  line  such  that  the  projected  distance  was  minimized.  Subsequently  the  station 
was  moved  along  the  line  until  its  distance  from  the  center  of  the  ridge  was  preserved.  The 
requirement  of  conservation  of  distance  from  the  ridge  was  linearly  relaxed  between  the 
center  of  the  basin  and  the  ridge  such  that  the  projection  was  completely  determined  by 
distance  from  the  ridge  near  the  ridge  and  completely  determined  by  minimum  projected 
distance  near  the  center  of  the  basin. 

In  total,  519  stations  from  the  32  sections  were  projected  onto  AR7E.  In  the  same  fash¬ 
ion  each  ADCP  measurement  was  projected  onto  the  AR7E  line.  The  only  difference  be¬ 
tween  these  procedures  was  that  the  original  depth  of  each  ADCP  measurement  was  solely 
determined  by  Etopo2.  Measurements  from  locations  requiring  a  projection  in  excess  of 
200  km  were  discarded.  As  most  of  the  hydrographic  stations  were  located  in  the  vicinity 
of  AR7E,  the  distance  of  projection  was  less  than  5  km  for  73%  of  the  stations  and  only 

^The  station’s  original  depth  was  taken  as  the  deepest  of  maximum  measurement  depth  plus  10  m  and 
Etopo2  interpolated  to  the  position  of  the  station 


77 


15%  of  the  stations  were  moved  more  than  25  km  from  their  original  locations  (Fig.  4-7). 
In  summary,  measurements  from  the  western  side  of  the  Irminger  Basin  were  projected 
along  isobaths  onto  the  AR7E  line,  while  in  the  eastern  half  of  the  basin  the  distance  from 
the  Reykjanes  Ridge  was  taken  into  account. 


Figure  4-7:  Bathymetry  profile  and  hydrographic  data  coverage  (in  percent).  The  black  triangles 
indicate  station  locations  along  AR7E.  The  inset  histogram  shows  the  distribution  of  distances  of 
projection. 

The  hydrographic  and  direct  velocity  data  from  individual  sections  were  then  interpo¬ 
lated  onto  a  standard  5  km  by  10  m  grid  using  a  Laplacian- spline  interpolator  (see  Pickart 
and  Smethie,  1998,  for  details  of  this  procedure).  The  data  coverage  along  the  standard  grid 
is  shown  in  Fig.  4-7.  The  western  half  of  the  Irminger  Basin  was  sampled  during  the  major¬ 
ity  of  the  occupations,  while  hydrographic  measurements  were  conducted  across  the  entire 
basin  from  the  Greenland  shelf  to  the  Reykjanes  Ridge  during  75%  of  the  occupations.  The 
collection  of  gridded  sections  will  be  used  throughout  the  rest  of  the  study. 

4.3.6  ADCP-referenced  geostrophic  velocities 

Absolute  geostrophic  velocities  along  the  AR7E  section  were  computed  two  ways:  using 
the  shipboard  ADCP  velocity  data  in  the  upper  part  of  the  water  column,  and  using  the 
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AVISO  gridded  surface  velocities.  Since  the  time  period  covered  by  AVISO  (1993-2007) 
is  much  greater  than  that  sampled  by  the  Nuka  Arctica  container  vessel  (1999-2002),  in 
order  to  meaningfully  investigate  interannual  variability  it  is  necessary  to  use  the  former 
product.  As  such,  the  ADCP-referenced  absolute  velocities  were  used  only  to  demonstrate 
the  accuracy  of  the  AVISO-referenced  velocities. 

For  the  Nuka  Arctica  time  period,  a  weighted  mean  ADCP  section  over  the  depth  inter¬ 
val  50  to  350  m  was  computed.  The  weights  are  inversely  proportional  to  the  distance  of 
projection  (projection  distances  shorter  than  25  km  were  given  equal  weight).  The  weight¬ 
ing  scheme  was  implemented  because  of  the  change  in  direction  of  the  current  flowing 
along  the  western  boundary  as  it  starts  to  curve  around  Cape  Farewell,  and  the  Greenland 
slope  is  also  the  only  region  where  the  weighted  and  unweighted  mean  vertical  ADCP 
sections  differed  significantly  (Fig.  4-8).  The  absolutely-referenced  geostrophic  flow  was 
computed  following  the  procedure  of  Pickart  et  al.  (2005).  In  particular,  each  of  the  eleven 
gridded  vertical  sections  of  temperature  and  salinity  for  the  time  period  1999-2002  was 
used  to  calculate  the  dynamic  height,  from  which  relative  geostrophic  velocities  were  com¬ 
puted.  The  individual  sections  were  then  re-mapped  onto  the  standard  grid  and  averaged 
to  create  a  mean  thermal  wind  section.  Finally,  this  section  was  referenced  by  match¬ 
ing  the  vertically  averaged  relative  and  ADCP  velocities  at  each  horizontal  grid  node  (the 
difference  between  a  vertically  averaged  match  and  a  least  squares  best  fit  was  minimal). 
Although  the  de-tided  ADCP  velocity  field  still  contains  ageostrophic  motion,  good  agree¬ 
ment  between  the  mean  ADCP  and  mean  absolute  geostrophic  velocities  (not  shown)  in¬ 
dicate  that  the  ageostrophic  components  were  to  a  great  extent  averaged  out  and  had  little 
influence  on  the  result.  Below  we  compare  the  near-surface  ADCP-referenced  geostrophic 
velocities  to  the  AVISO  product  for  the  4-year  period  of  overlap. 

4.3.7  AVISO-referenced  geostrophic  velocities 

The  mean  absolute  geostrophic  velocities  computed  from  the  subset  of  1 1  hydrographic  and 
the  47  direct  velocity  sections  provide  a  realistic  estimate  of  the  flow  across  the  AR7E  line 
during  the  1999-2002  period.  However,  to  examine  the  evolution  of  the  flow  throughout 
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Figure  4-8:  Comparison  of  ADCP  velocities  across  AR7E  for  the  period  1999-2002.  a)  Weighted 
mean  velocities,  b)  unweighted  mean  velocities,  and  c)  absolute  difference  between  a)  and  b). 


80 


Absolute  ADCP  velocity  difference  (cm/s) 


the  entire  sampled  period  (1991-2007),  additional  reference  velocity  information  is  needed. 
Global  sea  surface  height  has  continuously  been  measured  by  multiple  satellite  altimeters 
since  October  1992.  We  examined  the  merged  AVISO  geostrophic  surface  velocity  prod¬ 
uct  to  see  how  it  compares  against  our  independently  computed  ADCP-referenced  surface 
geostrophic  velocities.  The  mean  1999-2002  AVISO  surface  velocity  field  (for  the  months 
of  April  through  December,  to  coincide  with  the  period  of  hydrographic  and  direct  velocity 
measurements)  was  interpolated  onto  the  AR7E  line  and  the  velocity  component  normal  to 
AR7E  determined.  The  resulting  AVISO  velocity  curve  (green  line)  and  the  near-surface 
absolute  geostrophic  velocities  from  the  ADCP-referenced  field  (10-50  m  vertically  av¬ 
eraged,  blue  line)  are  displayed  in  Eigure  4-9a.  Although  the  AVISO  curve  is  smoother, 
the  agreement  between  the  AVISO-  and  ADCP-referenced  surface  geostrophic  velocities 
is  striking,  with  a  root-mean-square  (RMS)  deviation  of  less  than  3.0  cmj s.  Both  curves 
display  the  primary  circulation  features  of  the  Irminger  Sea:  the  Western  Boundary  Current 
(WBC,  approximately  0-200  km),  the  northward  branch  of  the  Irminger  Gyre  (IG,  approx¬ 
imately  200-350  km),  and  the  Irminger  Current  (IC,  approximately  500-700  km).  We  note 
that  the  AVISO  product  has  smoothed  out  the  two  surface  peaks  of  the  IC  that  were  pointed 
out  by  Knutsen  et  al.  (2005)  using  the  Nuka  Arctica  ADCP  data. 

Another  pair  of  surface  velocity  curves  computed  from  independent  mean  DOT  prod¬ 
ucts  were  evaluated  for  possible  use  as  reference  velocities  with  the  geostrophic  calculation 
(Pig.  4-9b).  These  would  need  to  be  coupled  to  the  AVISO  velocity  anomaly  fields  in  order 
to  provide  time- varying  reference  information.  The  GRACE-based  DOT  estimate  of  Tap- 
ley  et  al  (2003,  orange  line)  is  best  suited  for  applications  with  wavelengths  greater  than 
500  km,  and  does  not  properly  resolve  the  scales  of  interest  in  this  study.  Thompson  et  al 
(2009,  red  line)  used  terrestrial  and  altimeter-derived  gravity  data  to  add  high  wave  number 
detail  to  their  regional  DOT  of  the  northwest  Atlantic.  They  were  able  to  resolve  oceano¬ 
graphic  features  with  length  scales  of  order  100  km  such  as  the  Mann  Eddy,  and  do  also 
capture  the  primary  features  of  interest  in  the  Irminger  Sea  (the  WBC,  the  IG,  and  the  IC). 
However,  the  magnitudes  of  the  surface  flow  in  both  the  WBC  and  the  northward  branch  of 
the  IG  are  unrealistically  high,  and  the  peak  of  the  IC  is  shifted  eastward  by  about  100  km 
compared  to  our  ADCP-referenced  geostrophic  surface  velocities  (Pig.  4-9b).  A  follow-up 
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Figure  4-9:  Comparison  of  surface  velocities  across  the  AR7E  line,  a)  The  blue  line  is  a  trace 
of  the  mean  1999-2002  near-surface  absolute  geostrophic  velocities.  The  green  line  shows  the 
mean  1999-2002  AVISO  surface  velocities.  The  approximate  locations  of  the  primary  circulation 
features  (discussed  in  the  text)  are  indicated,  b)  The  additional  red  and  orange  lines  are  mean 
surface  velocities  computed  from  the  dynamic  ocean  topographies  (DOT)  of  Thompson  et  al.  (2009) 
and  Tapley  et  al.  (2003),  respectively. 


study  in  which  additional  high-quality  terrestrial  gravity  data  from  the  Irminger  Sea  area 
are  ineluded  will  be  eondueted  within  the  next  few  months  (K.  Thompson,  personal  com- 
munieation,  2009).  Despite  the  above  mentioned  shorteomings  (laek  of  aeouracy  for  flow 
on  the  shelf  and  on  short  seales),  the  AVISO  absolute  geostrophie  surfaee  veloeities  ealeu- 
lated  using  the  DOT  of  Rio  and  Hernandez  (2004)  are  obviously  the  best  choice  at  present 
for  refereneing  our  geostrophic  velocity  sections  over  the  full  time  period.  The  scales  of 
the  primary  features  of  interest  (the  WBC,  IG,  and  IC)  appear  to  be  properly  resolved,  and 
the  Greenland  shelf  area  will  not  be  ineluded  in  our  transport  estimates. 

In  order  to  investigate  the  time- varying  nature  of  the  flow  in  the  Irminger  Sea,  we  ehose 
to  eonsider  annual  averages,  or,  more  preeisely,  summer  averages  (May-Oet).  The  ratio¬ 
nale  was  to  reduee  the  effeets  of  mesoseale  aetivity  as  well  as  altimeter  data  uncertainty, 
and  the  majority  of  the  hydrographic  sections  were  oeeupied  during  these  months.  Accord¬ 
ingly,  a  single  composite  AVISO  veloeity  curve  was  eomputed  for  each  year  (May  through 
Oetober,  26-27  realizations)  for  the  period  1993-2007  (Fig.  4- 10a).  These  surfaee  velocity 
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curves  were  then  used  to  referenee  the  annually-averaged  thermal  wind  sections  calculated 
from  the  AR7E  oecupations.  AVISO  only  extends  baek  to  Oetober  1992,  whieh  is  not 
suffieient  to  provide  referenee  veloeity  information  for  the  1991  and  1992  hydrographie 
seetions.  However,  the  early  1990s  was  a  period  of  uniformly  strong  atmospherie  foreing, 
as  witnessed  by  the  eonsistently  high  NAO  index  (Fig.  4-4),  and,  perhaps  not  surprisingly, 
there  was  little  differenee  in  the  referenee  veloeity  eurves  from  1993  to  1996.  This  provided 
motivation  to  extend  the  timeseries  back  in  time  to  1991  via  linear  extrapolation  under  the 
assumption  that  the  eireulation  in  1991  and  1992  was  similar  in  eharaeter  to  the  sueeeed- 
ing  years  (Fig.  4-lOb)^.  This  is  admittedly  a  big  assumption,  and  the  resulting  transport 
estimates  for  1991  and  1992  should  be  eonsidered  as  less  eertain. 
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Figure  4-10:  Annually  (May-Oct)  averaged  AVISO  reference  velocity  curves,  a)  1993-2007.  b) 
1993-1996  (blue  lines)  and  1991-1992  extrapolated  surface  velocities  (red  lines). 


4.3.8  Absolute  geostrophic  velocity  errors 

It  is  important  to  assess  the  errors  assoeiated  with  the  ealeulation  of  absolute  geostrophie 
veloeities  refereneed  by  AVISO  surfaee  veloeities.  The  uncertainties  originate  both  from 

^Though  the  atmospheric  state  had  shifted  by  1996,  the  surface  circulation  had  not  yet  adjusted  to  make 
it  significantly  different  from  the  preceding  years,  and  it  was  included  as  the  most  recent  end  point  of  the 
extrapolation. 
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the  altimeter  and  the  hydrographic  data. 

One  important  difficulty  in  using  altimetry-derived  surface  velocities  as  reference  val¬ 
ues  for  full-depth  geostrophic  calculations  is  the  uncertainty  of  the  marine  geoid,  in  partic¬ 
ular  for  areas  in  the  vicinity  of  land  and  for  scales  shorter  than  about  500  km  (e.g.  Tapley 
et  al,  2003).  In  spite  of  this  we  found  very  good  correspondence  between  the  4-year  mean 
surface  velocities  derived  from  the  geostrophic  calculation  referenced  by  direct  velocity 
measurements  (ADCP)  and  the  AVISO  absolute  velocities  (Fig.  4-9).  This  indicates  that, 
averaged  over  a  sufficient  amount  of  time,  the  AVISO  surface  velocities  provide  a  reason¬ 
able  representation  of  the  directly  observed  surface  flow  across  the  AR7E  line.  From  this 
comparison  we  estimate  an  error  associated  with  the  geoid  based  on  the  RMS  difference 
between  the  ADCP-referenced  and  AVISO  geostrophic  surface  flows.  The  mean  RMS  dif¬ 
ference  is  3.0  cm/ s,  with  the  greatest  values  found  over  the  Greenland  shelf  and  slope  and 
smallest  values  in  the  region  of  the  Irminger  Gyre.  A  curve  of  RMS  differences  calculated 
from  a  running  mean  filter  with  a  window  of  200  km  provides  an  estimate  of  the  uncertainty 
of  the  geoid  at  each  nodal  point. 

Standard  errors  were  computed  for  each  of  the  summer  averaged  AVISO  reference  ve¬ 
locity  curves  (Fig.  4- 10a)  using  the  bootstrap  technique,  which  is  a  procedure  that  involves 
random  sampling  with  replacement  from  the  dataset  and  does  not  require  any  assumptions 
about  the  underlying  probability  distribution.  The  mean  standard  error  was  1.2  cm/ s,  and 
varied  between  approximately  0.8  cm/ s  in  the  WBC  and  1.6  cm/ s  in  the  IC.  The  largest 
standard  errors  were  generally  found  near  the  center  of  the  basin  and  in  the  IC.  We  assigned 
a  much  larger  error  to  the  extrapolated  1991  and  1992  reference  velocity  curves  estimated 
from  the  sum  of  the  range  and  mean  standard  error  of  the  1993-1996  curves  that  served  as 
the  basis  for  the  extrapolation. 

Density  gradients  cannot  be  directly  measured  below  the  deepest  common  level  of 
neighboring  hydrographic  stations  (i.e.  the  “bottom  triangle”).  The  extrapolation  resulting 
from  our  Laplacian- spline  gridding  technique  filled  these  bottom  triangles,  but  the  differ¬ 
ence  between  the  extrapolated  grid  points  and  the  real  ocean  is  unknown.  To  estimate  an 
error  associated  with  flow  through  the  bottom  triangles  we  first  determined  the  median  dis¬ 
tance  between  stations  (36  km)  and  the  median  slope  of  the  bathymetry  over  this  distance 
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(4.1  m/km)  to  obtain  a  “standard  bottom  triangle  area”.  In  the  vicinity  of  steep  bathymetry 
this  area  would  tend  to  be  larger.  But,  as  a  matter  of  routine,  station  spacing  is  typically  de¬ 
creased  to  compensate  for  the  increased  bottom  triangle  area  under  such  conditions,  thereby 
reducing  the  impact  of  steep  bathymetry  on  the  velocity  error.  The  resulting  standard  bot¬ 
tom  triangle  area  is  equivalent  to  a  uniform  layer  of  thickness  74  m  situated  directly  above 
the  bottom.  Combining  the  standard  deviation  of  the  absolute  geostrophic  velocities  within 
this  layer  and  the  area  of  a  standard  bottom  triangle,  we  obtain  a  transport  error  estimate  of 
0.18  5v  associated  with  the  data  gap.  Naturally  this  error  affects  only  transport  estimates  of 
currents  extending  to  the  bottom.  With  a  bottom  depth  of  2300  m  (the  average  depth  of  the 
Irminger  Basin)  and  a  median  distance  of  36  km  between  stations,  this  error  is  comparable 
to  a  barotropic  flow  of  0.22  cm/ s.  For  swift,  bottom-trapped  currents  such  as  the  DWBC, 
the  bottom  triangles  may  locally  represent  a  much  larger  uncertainty. 

An  error  will  also  result  from  the  projection  and  gridding.  Given  the  short  distance  over 
which  most  of  the  stations  were  projected  (85%  less  than  25  km)  we  do  not  expect  a  signif¬ 
icant  contribution  to  the  total  error  resulting  from  the  projection.  The  error  arising  from  the 
gridding  procedure  was  estimated  from  a  comparison  between  the  observed  data  and  the 
gridded  data  interpolated  to  the  locations  of  the  original  observations  (e.g.  Nikolopoulos 
et  al,  2009).  Below  250  m  the  median  RMS  difference  between  actual  observations  and 
gridded  maps  were  0.012°C  and  0.0021  for  temperature  and  salinity,  respectively.  The  tem¬ 
perature  RMS  difference  exceeds  the  observational  error  by  a  factor  of  6.  In  the  worst  case 
scenario,  where  one  station  would  consistently  be  biased  towards  colder,  saltier  conditions 
by  the  above  RMS  differences  compared  to  a  neighboring  station  separated  by  the  median 
distance  of  36  km,  the  geostrophic  velocity  error  would  reach  1  cm/ s  at  3000  m.  This  ex¬ 
treme  case  is  highly  unrealistic  due  to  the  random  component  of  the  gridding  error.  The 
greatest  RMS  difference  between  observed  and  gridded  data  occurred  on  the  Greenland 
shelf,  both  because  of  strong  gradients  in  temperature  and  salinity  as  well  as  large  differ¬ 
ences  in  the  station  spacing  between  the  various  cruises.  This  is  not  an  important  issue  for 
our  study,  since  the  Greenland  shelf  is  not  considered. 

Finally  there  is  a  statistical  error  associated  with  use  of  the  synoptic  density  structure  as 
a  proxy  for  a  summertime  mean  density.  In  the  Labrador  Sea  this  error  is  negligible  below 
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1500  m  (Han  and  Tang,  2001),  which  is  likely  the  case  for  the  Irminger  Sea  as  well.  This 
error  is  of  course  reduced  for  the  years  when  multiple  transects  are  available  (Fig.  4-4). 
We  have  not  attempted  to  quantify  this  statistical  error.  Instead  our  annual  mean  sections 
may  be  considered  synoptic  representations  for  each  given  year  smoothed  by  the  summer 
averaged  AVISO  reference  velocities. 

We  found  above  the  hydrographic  uncertainties  to  be  negligible  compared  to  the  errors 
associated  with  the  temporally  averaged  AVISO  absolute  surface  velocities  that  were  used 
for  referencing,  with  the  possible  exception  of  swift,  bottom-trapped  currents  near  steep 
bathymetry  (e.g.  the  DWBC)  where  the  error  arising  from  transport  through  the  bottom 
triangles  may  become  more  important.  The  total  uncertainty  of  the  reference  velocity  curve 
for  each  year  was  determined  as  the  root  of  the  sum  of  the  squares  of  the  two  dominant 
contributions  (geoid  and  standard  errors).  The  transport  errors  in  Section  4.6  are  estimated 
from  the  impact  of  the  reference  velocity  errors  on  the  geostrophic  calculations.  Velocities 
exceeding  ±  2  cm/5  may  in  general  be  considered  significant  for  any  pentadal  or  longer 
average. 


4.4  Mean  structure 

We  first  examine  the  overall  (1991-2007)  mean  hydrographic  and  absolute  velocity  fields, 
which  are  displayed  in  Figure  4-11.  Following  previous  studies  (e.g.  Holliday  et  ai,  2009) 
we  split  the  water  column  into  shallow,  intermediate,  and  deep  water  masses. 

Throughout  most  of  the  shallow  layer  (oe  <  21 .10  kg/ m^)  is  a  warm,  saline  Atlantic- 
origin  water  mass  called  Irminger  Sea  Water  that  is  mainly  supplied  by  the  Irminger  Current 
(IC).  The  IC  is  a  branch  of  the  North  Atlantic  Current,  and  enters  the  Irminger  Sea  as  a 
northward  flow  along  the  western  side  of  the  Reykjanes  Ridge.  In  the  mean  it  is  a  surface- 
intensified  flow  roughly  200  km  wide,  with  a  deep  extension  reaching  the  1500  m  isobath 
(Figure  4-11).  The  major  fraction  of  the  IC  (up  to  85%,  Pickart  et  al,  2005)  flows  around 
the  Irminger  Basin  (Fig.  4-2)  and  continues  south  along  the  Greenland  slope.  Here  a  sharp 
hydrographic  front  separates  the  Irminger  Sea  Water  from  the  cold,  fresh  Arctic-origin 
water  carried  by  the  East  Greenland  Current.  We  are  unable  to  distinguish  the  boundary 
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between  the  southward  flow  of  the  East  Greenland  Current  and  the  IC.  This  was  the  case  as 
well  in  studies  farther  to  the  south  (Fratantoni  and  Pickart,  2007).  We  will  hereafter  refer  to 
the  combined  southward  flow  along  the  Greenland  slope  as  the  Western  Boundary  Current 
and  the  northward  flow  along  the  Reykjanes  Ridge  as  the  Irminger  Current. 

The  intermediate  layer  (27.70  kg / <  Oq  <  27.80  kg/m^)  is  dominated  by  Labrador 
Sea  Water,  which  is  identifiable  by  a  mid-depth  minimum  of  salinity  and  potential  vortic- 
ity  (e.g.  Talley  and  McCartney,  1982).  The  low  potential  vorticity  signifies  weak  stratifica¬ 
tion,  which  may  be  traced  to  the  convective  origin  of  LSW  resulting  in  a  nearly  homoge¬ 
neous  water  mass.  The  most  pure  LSW  in  the  Irminger  Sea  is  found  in  the  western  half  of 
the  basin  (note  the  doming  and  deepening  27.70  and  27.80  kg/n?  isopycnals,  respectively). 
The  warmer  and  saltier  water  mass  occupying  the  same  depth  as  LSW  predominantly  in 
the  vicinity  of  the  Reykjanes  Ridge  is  called  Icelandic  Slope  Water  (van  Aken  and  de  Boer, 
1995;  Read,  2001). 

Two  distinct  water  masses  are  found  in  the  deep  layer  (Oe  >  27.80  kg/n?),  both  of 
which  contribute  to  the  DWBC  and  originate  from  the  dense  Nordic  seas  overflows.  Near 
the  bottom  the  cold  and  fresh  signature  of  the  Denmark  Strait  Overflow  Water  (DSOW)  is 
evident.  Directly  above  this  the  Iceland- Scotland  Overflow  Water  (ISOW)  is  identifiable 
by  the  deep  salinity  maximum.  This  water  mass  is  often  referred  to  as  North  East  Atlantic 
Deep  Water  (NEADW),  and  its  greater  salinity  results  largely  from  entrainment  of  saline 
Atlantic-origin  waters  primarily  in  the  vicinity  of  the  Iceland-Scotland  sills  (Fogelqvist 
et  al,  2003).  The  boundary  separating  ISOW  and  DSOW  is  usually  taken  to  be  the  27.88 
kgjn?  isopycnal  (e.g.  Dickson  and  Brown,  1994),  and  we  adopt  this  definition  as  well. 
Collectively  DSOW,  ISOW,  and  LSW  form  the  major  fraction  of  the  DWBC,  and  are  re¬ 
ferred  to  as  North  Atlantic  Deep  Water  (NADW),  a  water  mass  that  constitutes  a  significant 
fraction  of  the  World  Ocean  (Johnson,  2008)  and  plays  an  important  role  in  the  Meridional 
Overturning  Circulation. 

The  strongest  velocity  signature  in  the  mean  section  is  the  boundary  current  system. 
The  IC  flows  northward  on  the  western  side  of  the  Reykjanes  Ridge,  and  returns  as  a 
surface-intensified  flow  in  the  upper  part  of  the  WBC.  The  enhanced  southward  flow  at 
the  bottom  of  the  WBC  is  the  DSOW,  which  has  a  “tail”  extending  offshore  to  roughly  the 
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Figure  4-11:  Vertical  sections  of  a)  potential  temperature  (°C),  b)  salinity,  and  c)  absolute 
geostrophic  velocity  (positive  to  the  northeast,  cmjs)  averaged  over  the  period  1991-2007.  The 
grey  lines  are  contours  of  constant  potential  density,  and  the  black  and  green  lines  in  c)  are  the  zero 
velocity  contour  and  the  Irminger  Gyre  boundary.  The  black  triangles  indicate  station  locations. 
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3000  m  isobath.  This  feature  is  seen  in  synoptic  sections  as  well  (Pickart  et  al,  2005).  As 
previously  noted  (e.g.  Lherminier  et  al,  2009),  in  the  central  Irminger  Basin  there  is  little 
density  shear  through  most  of  the  water  column,  which  results  in  a  large  depth-independent 
component  to  the  flow.  The  northward  flow  of  the  Irminger  Gyre  (IG)  is  seen  clearly  in  the 
mean  section,  consistent  with  current  meter  data  collected  at  this  location  (Vage  et  al, 
2008).  We  defined  the  southward  limb  of  the  gyre  as  the  seaward  part  of  the  WBC  nec¬ 
essary  to  balance  the  northward  transport  at  each  density  level  in  the  water  column.  The 
boundary  of  the  IG  is  depicted  in  Figure  4-1  Ic  by  the  thick  green  line.  We  note  that  earlier 
studies  have  referred  to  the  Irminger  Gyre  as  the  cyclonic  circulation  across  the  full  width 
of  the  Irminger  Basin  (e.g.  Flatau  et  al,  2003;  Lherminier  et  al,  2007).  We  strictly  de¬ 
note  the  Irminger  Gyre  as  the  narrow  cyclonic  recirculation  in  the  western  Irminger  Basin. 
This  feature  was  first  mapped  by  Lavender  et  al.  (2000)  and  is  displayed  in  Figure  4-1. 
The  recirculating  component  of  the  WBC  forming  the  southward  limb  of  the  gyre  will  be 
removed  from  our  estimates  of  net  equatorward  WBC  transport. 

4.5  Temporal  variability 

Significant  changes  have  taken  place  across  the  Irminger  Basin  in  concert  with  the  chang¬ 
ing  atmospheric  conditions  represented  by  the  NAO  index  (Fig.  4-4).  To  document  these 
changes  we  first  focus  on  two  pentads,  the  early  1990s  (1991-95)  and  the  early  2000s 
(2001-05),  that  represent  distinctly  different  modes  of  the  atmospheric  circulation  as  rep¬ 
resented  by  the  NAO  index.  Following  this  we  relate  results  from  an  EOF  analysis  using 
hydrographic,  altimetric,  and  atmospheric  reanalysis  data  to  the  NAO  index  and  the  subpo¬ 
lar  gyre  index  of  Hakkinen  and  Rhines  (2004,  2009).  The  goal  is  to  understand  how  these 
changes  were  interrelated. 

4.5.1  Contrasting  the  early  1990s  to  the  early  2000s 

The  early  1990s  (1991-95)  was  a  period  of  consistently  high  NAO  index.  It  was  charac¬ 
terized  by  stormy  and  severe  winters,  and  the  formation  of  a  large  volume  of  LSW  in  the 
Labrador  Sea  (Lazier  et  al,  2002;  Yashayaev,  2007).  The  resulting  LSW  class,  the  coldest 
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and  densest  on  record,  has  been  referred  to  as  lower  LSW  (ILSW)  or  LSW1994.  Following 
this,  the  NAO  index  dropped  precipitously  to  a  strongly  negative  value  in  1996,  and  the 
formation  of  ILSW  ceased.  In  the  early  2000s  (2001-05)  the  NAO  index  fluctuated  be¬ 
tween  weakly  positive  and  weakly  negative.  The  atmospheric  forcing  at  the  beginning  of 
the  period  (winter  1999-00)  was  sufficient  to  cause  renewed  formation  of  LSW,  although  a 
less  dense  class  of  LSW  was  formed,  called  upper  LSW  (uLSW)  or  LSW2000  (Yashayaev, 
2007;  Yashayaev  et  al,  2007a). 

To  illustrate  the  changes  that  took  place  between  these  vastly  different  periods  of  atmo¬ 
spheric  forcing,  we  display  in  Figures  4-12  and  4-13  mean  vertical  sections  for  the  early 
1990s  and  early  2000s.  The  most  dramatic  change  in  the  hydrography  of  the  Irminger  Sea 
occurred  in  the  intermediate  layer  (Figs.  4-12  and  4-13).  A  large  pool  of  cold,  fresh,  and 
weakly  stratified  ILSW  filled  most  of  the  western  and  central  Irminger  Basin  between  the 
27.75  and  27.80  kg/n?  isopycnals  in  the  early  1990s.  By  the  early  2000s  the  ILSW  dome 
had  all  but  collapsed,  and  warmer,  saltier,  and  more  stratified  waters  permeated  the  basin. 
As  the  formation  of  ILSW  ceased  in  the  mid-1990s,  the  water  mass  became  completely 
isolated  from  the  surface.  Mixing  with  warmer  and  saltier  waters  along  the  margins  eroded 
the  salinity  minimum  and  led  to  restratification  and  flattening  of  the  isopycnals  (e.g.  Lazier 
et  al.,  2002).  In  the  early  2000s  there  was  still  a  trace  of  ILSW  identified  by  the  deep  salin¬ 
ity  minimum  above  the  27.80  kg/n?  isopycnal.  The  shallower  salinity  minimum  above  the 
27.75  kg/n?  isopycnal  is  associated  with  uLSW.  These  large-scale  changes  in  the  lower 
and  upper  LSW  volumes  throughout  the  subpolar  gyre  are  readily  detectable  in  chloroflu¬ 
orocarbon  inventories  (e.g.  Kieke  et  al,  2007).  The  local  salinity  maximum  separating 
the  two  LSW  modes  is  referred  to  as  Icelandic  Slope  Water  (Holliday  et  al,  2009).  The 
weak  stratification  in  the  vicinity  of  the  Reykjanes  Ridge  is  also  a  signature  of  this  water 
mass  (Read,  2001),  and  is  not  related  to  LSW.  The  importance  of  local  formation  of  LSW 
in  the  Irminger  Sea  versus  import  from  the  Labrador  Sea  will  be  discussed  in  Section  4.7. 
Warming  and  salinification  took  place  in  the  deep  layer  as  well,  and  has  been  quantified  by 
previous  authors  (Sarafanov  et  al,  2007;  Yashayaev  et  al,  2007a,b). 

The  declining  buoyancy  forcing  and  collapse  of  the  ILSW  dome  during  the  1990s  re¬ 
sulted  in  a  rising  subpolar  sea  surface  height  and  a  spin-down  of  the  subpolar  gyre  (Hakkinen 
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black  triangles  indicate  station  locations. 
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Figure  4-13:  As  in  Figure  4-12,  but  for  a,  b)  potential  vorticity  {{ms)  *  *  10  *^)  and  c,  d)  absolute 
geostrophic  velocity  {cm/s). 
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and  Rhines,  2004).  This  is  reflected  in  the  changing  structure  of  the  boundary  currents 
along  the  margins  of  the  Irminger  Sea.  The  shallow  part  of  the  IC  weakened  between  the 
early  1990s  and  the  early  2000s.  The  most  notable  change  in  the  IC,  however,  is  in  its 
deep  structure.  During  the  latter  pentad  a  second  deep  extension  of  the  current  developed, 
centered  on  the  2500  m  isobath  (there  is  a  weak  signature  of  this  feature  during  the  first 
pentad).  The  influence  of  this  varying  IC  structure  on  the  transport  of  the  current  will  be 
quantified  in  Section  4.6,  and  possible  reasons  for  the  structural  transformation  will  be  dis¬ 
cussed  in  Section  4.8.  The  magnitude  of  the  WBC  also  declined  between  the  early  1990s 
and  the  early  2000s,  in  particular  in  the  shallow  layer,  consistent  with  the  notion  of  a  de¬ 
clining  subpolar  gyre.  Interestingly  there  is  a  compensating  effect  in  that  the  WBC  became 
wider  in  the  latter  period.  Timeseries  of  transport  in  the  WBC  will  also  be  examined  in  Sec¬ 
tion  4.6.  An  unexpected  change  that  occurred  between  the  two  periods  is  that  the  strength 
of  the  IG  (located  between  approximately  150  and  350  km  along  AR7E)  increased.  This 
is  surprising  because  both  the  wind  stress  and  buoyancy  components  of  the  atmospheric 
forcing  were  significantly  reduced  over  this  time  period  (Fig.  4-14),  and  local  wind  stress 
curl  and  air-sea  buoyancy  exchange  are  thought  to  be  the  main  forcing  mechanisms  of  the 
IG  (Spall  and  Pickart,  2003;  Pickart  et  al,  2003a). 

Examination  of  the  AVISO  DOT  field  (Fig.  4-15)  provides  an  indication  of  what  might 
have  taken  place.  In  agreement  with  the  observations  of  Hakkinen  and  Rhines  (2004)  an 
increase  in  SSH  is  evident  across  the  Irminger  Sea  (Fig.  4-I5a).  More  detailed  scrutiny 
reveals,  however,  that  the  increase  was  not  uniform.  The  extent  of  sea  level  rebound  was 
greater  in  the  eastern  half  of  the  Irminger  Basin  than  in  the  western.  As  a  result,  the  SSH 
gradient  in  the  central  Irminger  Sea  steepened.  This  effect  is  particularly  evident  when 
examining  the  DOT  along  AR7E  in  Figure  4-16.  Consequently,  the  northward  flow  of 
the  IG,  which  is  largely  barotropic,  increased.  This  is  consistent  as  well  with  the  AVISO 
surface  geostrophic  velocity  field  (Fig.  4- 15c,  d),  where  enhanced  northward  velocities 
across  AR7E  in  the  region  of  the  IG  is  evident  in  2001-05  relative  to  1992-95.  The  velocity 
difference  field  (not  shown)  shows  increased  velocities  in  most  of  the  interior  Irminger  Sea 
north  of  AR7E.  Increased  surface  flow  is  also  evident  along  parts  of  the  IC  (although  not 
at  the  AR7E  line  as  discussed  above),  which  is  surprising  given  the  general  decline  of  the 
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Figure  4-14:  Maps  of  NARR  winter  a,  b)  wind  stress  curl  (N/m^  *  10^^)  and  wind  vectors  (m/s) 
and  c,  d)  total  turbulent  heat  flux  (W/m^)  averaged  over  the  periods  a,  c)  1991-1995  and  b,  d)  2001- 
2005.  The  black  line  and  numbers  indicate  the  AR7E  line  and  distance  (in  km),  depth  is  contoured 
with  500  m  increments  starting  at  500  m. 
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subpolar  gyre  circulation.  The  WBC  surface  flow  declined  as  expected  from  Hakkinen  and 
Rhines  (2004).  The  apparent  paradox  of  a  strengthening  IG  circulation  under  weakening 
atmospheric  forcing  will  be  discussed  further  in  Section  4.8. 

4.5.2  Decadal  variability 

Hakkinen  and  Rhines  (2004)  used  empirical  orthogonal  functions  (EOFs)^  to  study  the 
North  Atlantic  sea  level  anomalies.  Interestingly,  the  subpolar  gyre  center  of  action  in 
the  spatial  pattern  of  their  leading  EOF  mode  is  situated  in  the  Irminger  Sea.  Hakkinen  and 
Rhines  (2004)  speculate  that  this  may  be  a  reflection  of  changes  in  local  convective  activity. 
Too  much  emphasis  should  perhaps  not  be  attached  to  the  exact  location  of  the  center  of 
action  considering  the  low  amount  of  variance  in  the  altimeter  data  set  explained  by  their 
leading  mode  (11%)  and  the  southern  Irminger  Sea’s  close  proximity  to  the  geographical 
center  of  the  subpolar  gyre. 

We  investigated  the  decadal  variability  with  a  more  local  focus,  and  performed  an  EOF 
analysis  using  a  wider  variety  of  oceanic  and  atmospheric  data  from  the  Irminger  Sea  re¬ 
gion.  In  particular,  we  computed  EOFs  using  the  AR7E  hydrographic  data,  altimetric  data, 
and  NARR  atmospheric  data.  We  performed  both  individual  and  coupled  EOFs,  including 
one  calculation  with  all  of  the  fields.  The  leading  modes  of  all  of  the  calculations  are  sim¬ 
ilar  in  character  and  represent  the  major  evolution  that  took  place  between  the  early  1990s 
and  the  early  2000s.  We  present  only  the  PC  timeseries  and  not  the  spatial  patterns,  which 
are  better  represented  by  the  pentadal  averages  presented  above  (Figures  4-12  to  4-15). 

For  the  hydrographic  EOF  we  used  as  input  the  annually-averaged  potential  tempera¬ 
ture,  salinity,  and  potential  vorticity  from  the  intermediate  layer  in  the  IG  region  to  empha¬ 
size  the  FSW.  The  resulting  PC  timeseries  explains  38%  of  the  variance  and  is  expressed 
by  the  black  crosses  in  Figure  4- 17a.  The  other  timeseries  in  the  figure  are  the  NAO  in¬ 
dex  (blue  line)  and  the  subpolar  gyre  index  (the  PC  timeseries  from  Hakkinen  and  Rhines 
(2004,  2009),  red  line).  We  note  immediately  the  high  level  of  agreement  between  the 

^Performing  an  EOF  analysis  of  a  data  set  implies  decomposing  it  into  basis  functions  orthogonal  to  each 
other,  chosen  such  that  the  variance  accounted  for  by  each  is  maximized.  Timeseries  and  spatial  patterns  for 
each  basis  function  are  both  produced  by  the  analysis. 
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Figure  4-15:  As  Figure  4-14,  but  for  AVISO  summer  absolute  a,  b)  dynamic  ocean  topography 
(cm)  and  c,  d)  surface  geostrophic  velocity  {cm/s).  The  red  segment  along  AR7E  indicates  the 
width  of  the  IG,  shown  also  in  Fig  3-15. 
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Figure  4-16:  AVISO  DOT  (cm)  along  AR7E  for  the  early  1990s  and  the  early  2000s.  The  width  of 
the  IG  is  indicated  by  the  red  segment. 

PC  timeseries  and  the  gyre  index.  The  correlation  between  the  two  timeseries  is  0.95.  It 
is  rather  remarkable  that  the  evolution  of  the  gyre  index,  which  represents  the  large-scale 
circulation  in  the  North  Atlantic,  so  closely  reflects  the  variability  of  the  Irminger  Sea  hy¬ 
drography  as  seen  in  the  AR7E  sections.  The  precise  timing  of  the  mid-1990s  maximum  is 
the  only  discrepancy  between  the  timeseries.  The  gyre  index  peaked  in  early  1994,  while 
the  PC  timeseries  peaked  around  1995.  The  difference  between  the  two  years,  on  the  other 
hand,  is  probably  insignificant,  and  the  missing  data  from  1993  would  likely  have  had  an 
effect  on  the  hydrographic  EOF  analysis.  With  the  exception  of  the  absolute  geostrophic 
velocity  field,  which  as  discussed  earlier  behaved  contrary  to  expectations  in  terms  of  the 
IG  circulation,  various  other  combinations  of  hydrographic  data  produce  similar  results. 
The  highest  correlation  with  the  gyre  index  (0.97)  was  obtained  when  computing  the  EOF 
from  potential  vorticity  alone,  and  the  greatest  amount  of  variance  explained  (53%)  was 
achieved  when  considering  potential  temperature  alone.  Expanding  the  domain  from  the 
intermediate  layer  in  the  IG  to  the  full  width  and  depth  of  the  Irminger  Sea  led  to  a  marginal 
increase  in  coherence  with  the  gyre  index,  but  a  significant  decrease  in  amount  of  variance 
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accounted  for.  Since  only  75%  of  the  sections  extended  across  the  entire  Irminger  Basin, 
fewer  data  points  were  available  from  the  full  domain. 

Next  we  considered  an  altimetric  EOF  using  the  summer  mean  dynamic  ocean  topog¬ 
raphy  and  surface  geostrophic  velocity  fields.  The  resulting  PC  timeseries  (Fig.  4- 17b) 
explains  28%  of  the  variance.  Although  the  EOF  is  computed  from  a  limited  domain  (the 
Irminger  Sea  deeper  than  1000  m,  north  of  57°N,  and  between  45°W  and  the  Reykjanes 
Ridge),  the  input  data  express  essentially  the  same  signal  as  the  sea  level  anomalies  used 
by  Hakkinen  and  Rhines  (2004,  2009),  and  the  high  correlation  (0.95)  between  the  PC  time- 
series  and  the  gyre  index  is  not  surprising.  The  coherence  is  greater  still  when  considering 
only  the  velocity  field,  while  most  variance  is  explained  by  the  first  DOT  EOF  mode. 

The  leading  atmospheric  EOF  mode  of  winter  mean  turbulent  heat  fluxes  and  wind 
stress  curl  from  the  NARR  atmospheric  reanalysis  project  (same  domain  as  above)  ex¬ 
plains  29%  of  the  variance  in  the  data  set.  Inspection  of  the  PC  timeseries  (Fig.  4-I7c) 
shows  a  high  degree  of  coherence  with  the  NAO  index,  which  indicates  that  the  NAO  in¬ 
dex  is  useful  as  a  representation  of  the  wintertime  atmospheric  forcing  in  this  region.  The 
strongest  signal  is  found  directly  east  of  Cape  Farewell  (Fig.  4-14),  which  is  a  testament 
to  the  importance  of  the  tip  jet  on  this  region  (e.g.  Pickart  et  al,  2003a).  As  noted  above 
there  is  also  a  positive  correlation  between  the  NAO  index  and  the  number  of  tip  jet  events 
per  winter.  The  region  influenced  by  the  atmospheric  variability  represented  by  the  NAO 
index  extends  well  beyond  the  Irminger  Sea.  Though  the  NAO  index  displays  a  substantial 
amount  of  short-term  fluctuations,  the  spin-down  and  contraction  of  the  subpolar  gyre  is 
likely  a  reflection  of  a  long-term  decline  in  the  NAO  index  from  its  maximum  in  the  early 
1990s  (see  for  example  Figure  Id  in  Sarafanov,  2009).  This  is  not  very  evident  from  the 
temporally  limited  PC  timeseries  in  Figure  4-17,  but  it  should  not  be  surprising  that  the 
ocean  and  the  atmosphere  operate  on  different  time  scales. 

4.6  Transports 

Transports  of  the  three  primary  features  of  interest,  the  Irminger  Gyre,  the  Irminger  Cur¬ 
rent,  and  the  Western  Boundary  Current,  have  been  estimated  using  the  annually  averaged 
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Figure  4-17:  Principal  components  timeseries  of  the  first  EOF  mode  computed  from  a)  annual 
mean  hydrographic  data  (potential  temperature,  salinity,  and  potential  vorticity),  b)  summer  mean 
altimeter  data  (dynamic  ocean  topography  and  surface  geostrophic  velocity),  and  c)  winter  mean 
atmospheric  reanalysis  data  (turbulent  heat  fluxes  and  wind  stress  curl).  The  blue  and  red  timeseries 
are  the  NAO  index  and  the  gyre  index  of  Hakkinen  and  Rhines  (2004,  2009),  respectively. 
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absolute  geostrophic  velocity  sections  described  above.  In  the  following  figures  (Figs.  4- 
18  to  4-21)  the  black  error  bars  represent  the  full  uncertainty.  This  is  a  combination  of  the 
geoid  error  estimate  and  the  AVISO  reference  velocity  standard  error,  which  are  the  two 
dominant  error  terms  (Sec.  4.3.8).  The  geoid  error  is  constant  in  time.  ^  The  blue  error  bars 
represent  only  the  AVISO  standard  error,  which  does  vary  from  year  to  year,  and  hence  is 
the  only  dominant  error  term  that  will  directly  impact  the  trend.  The  reader  should  keep  in 
mind  the  extrapolation  used  to  produce  absolute  geostrophic  velocity  sections  for  1991  and 
1992,  and  hence  the  increased  uncertainty  of  the  resulting  transports. 

4.6.1  The  Irminger  Gyre 

The  Irminger  Gyre  was  typically  located  between  150  and  350  km  along  the  AR7E  line 
(Fig.  4-1  Ic).  Although  it  was  present  every  summer,  there  was  considerable  interannual 
variability  in  terms  of  its  width  and  magnitude.  The  transport  was  estimated  from  the 
northeasterly  return  flow  to  the  east  of  the  center  of  the  gyre  (assuming  an  equal  return 
flow  in  the  WBC).  The  resulting  transport  timeseries  of  the  Irminger  Gyre  is  presented  in 
Figure  4-18.  As  noted  above  in  Section  4.5  the  magnitude  of  the  gyre  circulation  increased 
through  the  observational  period,  which  is  confirmed  by  the  trend  of  increasing  transport  at 
a  rate  of  4.3  Sv  per  decade.  The  significance  of  the  trend  was  found  to  be  greater  than  99% 
using  the  reverse  arrangements  test  of  Bendat  and  Piersol  (1986).^  The  mean  Irminger  Gyre 
transport  through  the  full  17-year  period  was  6.8  ±  1.8  5v.  An  alternative  way  to  compute 
the  mean  transport  of  the  IG  is  directly  from  the  temporally  averaged  absolute  geostrophic 
velocity  (Fig.  4-1  Ic)  instead  of  from  the  transport  timeseries.  The  resulting  mean  transport 
estimate  is  lower,  only  4.0  ±1.9  Sv.  This  discrepancy  exists  because  the  location  of  the 
gyre  changes  from  year  to  year  as  well  as  its  strength.  Hence  the  average  of  the  individual 
timeseries  values  can  be  thought  of  as  a  “stream  coordinate”  mean. 

^The  geoid  error  will  appear  to  have  temporal  variability  as  the  width  of  the  feature  under  consideration 
was  not  constant;  for  the  same  transport  a  broad  and  sluggish  current  will  have  a  greater  uncertainty  than  a 
narrow  and  swift  current. 

reversal  in  a  sequence  of  observations  of  the  variable  x  occurs  when  x,  >  Xj  for  i  <  j.  The  reverse 
arrangements  test  is  used  to  evaluate  the  significance  of  a  trend  based  on  how  many  such  reversals  take  place 
relative  to  the  number  expected  from  pure  chance  with  no  trend. 
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Figure  4-18:  Irminger  Gyre  transport  in  Sv.  The  black  error  bars  indicate  the  total  error,  while  the 
blue  error  bars  show  the  AVISO  standard  error  only.  The  red  line  is  the  least  squares  best  fit  through 
the  data. 


4.6.2  The  Irminger  Current 

The  Irminger  Current  is  steered  by  the  topography  of  the  Reykjanes  Ridge,  and  was  typi¬ 
cally  found  between  500  and  700  km  along  AR7E  (Fig.  4-1  Ic).  The  full-depth  IC  transport 
timeseries  is  displayed  in  Figure  4- 19a.  Surprisingly,  we  found  a  trend  of  increasing  IC 
transport,  which  is  contrary  to  the  notion  of  a  declining  subpolar  gyre  circulation.  There  is, 
however,  considerable  interannual  variability,  and  the  3.4  Sv  per  decade  trend  is  significant 
only  at  the  92%  level.  Removal  of  the  high-uncertainty  transport  estimates  from  1991  and 
1992  referenced  by  extrapolated  surface  velocities  does  little  to  improve  the  significance  of 
the  trend.  Closer  inspection  reveals  that  the  positive  trend  is  caused  by  a  transport  increase 
at  depth.  In  fact,  the  transport  in  the  upper  500  m  decreased  into  the  2000s  (Fig.  4- 19b), 
consistent  with  the  conclusions  of  Hakkinen  and  Rhines  (2004,  2009).  ^  The  deep  trans¬ 
port  increase  in  the  IC  (Fig.  4- 19c)  is  largely  the  result  of  the  emergence  of  the  second  deep 
extension  of  the  IC  noted  in  Section  4.5.  The  changing  structure  of  the  deep  IC  will  be 
discussed  more  extensively  in  Section  4.8,  including  possible  causes. 

The  mean  of  the  full-depth  IC  transport  timeseries  was  20  ±  3  Sv.  For  comparison  with 

®Both  the  weakening  trend  of  upper-layer  transport  and  the  strengthening  trend  of  the  deep  transport  were 
significant  at  better  than  the  95%  confidence  level. 
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Figure  4-19:  As  Figure  4-18,  but  for  the  Irminger  Current,  a)  Full  1C  transport,  b)  1C  transport 
above  500  m,  and  c)  1C  transport  below  1000  m.  Note  the  different  vertical  scale  in  b)  and  c). 


102 


previous  estimates  we  also  computed  the  transport  of  waters  less  dense  than  ae  =  27.76 
kg/rn^  whose  average  value  was  12  ±  2  5v.  The  corresponding  synoptic  transport  reported 
by  Pickart  et  al.  (2005)^  for  their  August  2001  occupation  was  13.7  Sv.  While  this  is  a 
bit  larger  than  our  average  value,  it  compares  well  with  our  2001  estimate  of  14  ±  6  5v. 
IC  transports  were  also  provided  by  Bacon  (1997):  9  Sv  (02  <  36.873  kgjn?  in  1991) 
and  Lherminier  et  al.  (2009):  1 1  ±  2  5v  (02  <  36.94  kg/rn^  in  both  2002  and  2004). 

The  IC  transport  computed  directly  from  the  temporally  averaged  vertical  section  of  ab¬ 
solute  geostrophic  velocity  (Fig.  4-1  Ic)  is  14  ±  3  Sv.  As  was  the  case  for  the  IG  transports, 
this  value  is  reduced  compared  to  the  temporal  average  of  the  discrete  estimates  presented 
above.  In  the  IC  case,  however,  meandering  of  the  current  is  not  the  cause  of  this  reduction; 
rather,  it  is  due  primarily  to  the  significant  changes  in  the  current’s  deep  structure,  where 
pentadal  flow  reversals  tend  to  reduce  the  deep  northward  transport  (compare  Figures  4-13 
c  and  d). 

4.6.3  The  Western  Boundary  Current  system 

We  consider  the  Western  Boundary  Current  (WBC)  to  be  the  southward  flow  bounded  by 
the  Greenland  shelfbreak  to  the  west  (30  km  along  AR7E)  and  the  deepest  part  of  the 
Irminger  Gyre  in  the  east  (typically  located  around  280  km,  see  Fig  3-10).  The  bottom  of 
the  Irminger  Gyre  (sometimes  it  extended  all  the  way  to  the  bottom)  is  close  to  the  deepest 
point  of  the  Irminger  Basin,  and  thus  our  definition  of  the  WBC  is  in  general  agreement  with 
that  of  Holliday  et  al.  (2009),  although  they  also  included  the  outer  part  of  the  Greenland 
shelf  in  their  calculations.  In  1995  the  section  did  not  extend  farther  onshore  than  to  50 
km.  To  compensate  for  the  missing  part  of  the  WBC  transport  that  year,  the  mean  transport 
between  30  and  50  km  (-3.2  Sv)  was  added  to  the  1995  estimate,  and  the  error  was  increased 
by  the  standard  deviation  of  the  mean  transport  within  the  given  interval  (1.9  Sv).  As  the 
DWBC  does  not  extend  inshore  of  50  km,  it  is  not  affected  by  the  modification  to  the  1995 
velocity  section. 

Traditionally,  all  of  the  southward  flow  on  the  western  side  of  the  Irminger  Basin  is 

^Pickart  et  al.  (2005)  also  provides  a  list  of  previous,  unreferenced  IC  transport  estimates,  none  of  which 
exceed  10  Sv. 
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taken  into  account  in  estimates  of  WBC  transport.  For  comparison  with  previous  results 
we  have  computed  WBC  transports  in  the  same  fashion,  hereafter  called  southward  WBC 
transport.  We  will,  however,  emphasize  timeseries  of  WBC  transport  with  the  IG  return 
flow  removed,  which  is  referred  to  as  the  net  WBC  transport.  The  net  WBC  transport 
is  the  throughput  of  the  boundary  current  system,  hence  the  more  appropriate  quantity  to 
consider. 

The  net  surface  to  bottom  WBC  average  transport  (Fig.  4-20)  was  -32  ±  4  5v  (keep  in 
mind  that  southward  transports  are  negative).  There  is  a  slight  indication  of  weakening  of 
the  WBC,  but  the  trend  is  not  significant.  The  range  of  interannual  variability  was  15  Sv. 
When  considering  the  southward  WBC  transport,  the  average  transport  is  larger,  -39  ±  4 
Sv.  Interestingly,  there  was  a  significant  trend  of  increasing  southward  WBC  transport  of 
-3.0  5v  per  decade.  However,  this  is  predominantly  because  the  Irminger  Gyre  spun  up  over 
this  time  period  (by  4.3  Sv  per  decade,  see  above).  Hence  there  was  little  or  no  trend  in  net 
transport.  In  the  case  of  the  WBC  there  was  good  agreement  between  the  mean  transport 
computed  from  the  transport  timeseries  and  from  the  temporally  averaged  vertical  section 
of  absolute  geostrophic  velocity  (-37  ±  3  Sv,  Fig.  4-1  Ic).  This  was  expected  given  the  lack 
of  WBC  spatial  variability  and  changing  structure,  unlike  the  IG  and  IC  cases  above. 
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Figure  4-20:  As  Figure  4-18,  but  showing  net  Western  Boundary  Current  transports. 


Our  southward  WBC  transport  (-39  ±  4  Sv)  compares  well  with  a  recent  estimate 
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by  Holliday  et  al.  (2009),  who  reported  a  transport  of  -40.5  ±  8.1  5v  in  August-September 
2005.  Earlier  surface  to  bottom  estimates  from  the  southern  Irminger  Sea  and  Cape  Farewell 
region  include:  -33.5  Sv  (2-month  current  meter  records  from  1978,  Clarke,  1984),  -32.2 
Sv  (ADCP-referenced  synoptic  section  from  1991,  Bersch,  1995),  -26.5  Sv  (inverse  model 
constrained  by  ADCP  data  from  1991,  Bacon,  1997),  and  -39.0,  -31.6,  and  -34.9  Sv  (also 
an  inverse  model  constrained  by  ADCP  data,  from  1997,  2002,  and  2004,  respectively, 
Lherminier  et  al,  2007,  2009).  Our  result  is  in  the  upper  range  of  the  previous  transport 
estimates. 

4.6.4  The  Deep  Western  Boundary  Current 

Given  its  importance  to  the  MOC,  we  have  considered  the  DWBC  component  (oe  >  27.80 
kg/m^)  of  the  WBC,  which  is  composed  of  the  DSOW  and  ISOW.  To  compensate  for  the 
bottom  triangle  error  (see  Section  4.3.8),  we  have  added  0.3  and  1.4  Sv  to  the  total  errors  of 
the  ISOW  and  DSOW  layers,  respectively  (estimated  from  the  average  0.18  Sv  per  36  km 
bottom  triangle  error  and  the  mean  distance  over  which  each  layer  is  in  contact  with  the 
bottom).  As  in  the  case  of  the  WBC  we  emphasize  net  transports  (although  for  the  DSOW 
layer,  into  which  the  IG  only  sporadically  reached,  the  difference  between  southward  and 
net  transports  is  minimal). 

The  net  DWBC  average  transport  (Fig.  4-21)  was  -9.6  ±1.3  Sv.  This  can  be  split  into 
-7.3  ±  1.0  Sv  and  -2.4  ±  0.4  Sv  for  the  ISOW  and  DSOW  components,  respectively.  Dur¬ 
ing  the  17-year  observational  period  the  net  DWBC  transport  increased  at  a  rate  of  -1.4 
Sv  per  decade  (the  trend  was  significant  at  a  level  higher  than  95%).  Using  southward 
DWBC  transport  estimates  relative  to  a  level  of  no  motion  at  1000  m,  Sarafanov  et  al. 
(2009a)  found  a  similar  significant  trend  of  increasing  baroclinic  transport  over  the  same 
time  period.  Combining  altimetry-derived  surface  velocities  and  baroclinic  velocities  from 
hydrographic  data  Sarafanov  et  al.  (2009b)  determined  that  the  absolute  DWBC  transport 
increased  as  well  between  the  mid-1990s  and  the  mid-2000s,  consistent  with  our  results. 
A  trend  of  increasing  DWBC  transport  in  the  Fabrador  Sea  from  the  late  1990s  to  the 
early  2000s  was  reported  by  Dangler  et  al.  (2006).  Extending  their  timeseries  using  equiv- 
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alently  computed  baroclinic  transports  from  Bacon  (1998),  a  robust  negative  correlation 
between  DWBC  transport  anomalies  and  LSW  thickness  anomalies  was  detected,  such  that 
decreasing  LSW  thickness  is  associated  with  enhanced  DWBC  transport.  This  relationship 
was  found  farther  south  (west  of  the  Grand  Banks  of  Newfoundland)  as  well  by  Pickart 
et  al.  (1999).  Sarafanov  et  al.  (2009a,b)  hypothesize  that  complex  and  coherent  oceanic 
responses  to  the  NAO-modulated  atmospheric  forcing  that  influences  (among  other  things) 
the  convective  activity  and  subpolar  gyre  circulation  may  regulate  the  relationship  between 
LSW  thickness  and  DWBC  transport. 


Figure  4-2 1 :  Net  Deep  Western  Boundary  Current  transport  in  Sv.  The  black  error  bars  indicate 
the  total  error  (including  estimated  bottom  triangle  error),  and  the  red  line  is  the  least  squares  best 
fit  through  the  data. 


For  the  full  DWBC  we  obtain  a  mean  southward  transport  of  -11.4  ±  1.2  5v.  This 
estimate  is  well  within  the  range  of  previous  synoptic  absolute  DWBC  transports  of  -7.3 
to  -13.3  Sv  for  the  southern  Irminger  Sea  and  Cape  Farewell  region  (e.g.  Clarke,  1984; 
Dickson  and  Brown,  1994;  Bersch,  1995;  Lherminier  et  al.,  2007,  2009;  Holliday  et  al., 
2009). 
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4.7  Convection 


The  traditional  hypothesis  describing  the  LSW  formation  process,  proposed  by  McCartney 
and  Talley  (1982),  argues  that  LSW  is  the  end-product  of  “subpolar  mode  waters”  that  are 
continuously  densified  via  air-sea  interaction  within  the  cyclonic  flow  around  the  subpo¬ 
lar  gyre.  This  culminates  with  deep  convection  in  the  Labrador  Basin,  forming  the  weakly 
stratified  LSW.  In  contrast,  a  growing  body  of  evidence  shows  that  deep  convection  can  also 
take  place  in  the  Irminger  Sea  during  particularly  severe  winters,  leading  to  a  second  loca¬ 
tion  of  LSW  formation  (e.g.  Pickart  et  al,  2003a,b;  Bacon  et  al,  2003;  Falina  et  al,  2007; 
Vage  et  al,  2008,  2009a).  These  studies  have  underscored  the  importance  of  small-scale 
oceanic  (the  Irminger  Gyre)  and  atmospheric  (the  westerly  Greenland  tip  jet)  features  that 
conspire  to  provide  conditions  favorable  for  deep  convection  in  the  southwestern  Irminger 
Sea,  and  suggest  that  high  spatial  resolution  may  be  necessary  to  understand  the  convective 
process  in  this  region. 

Despite  this  body  of  evidence,  recent  studies  argue  that  the  LSW  signature  in  the 
Irminger  Basin  is  principally  due  to  advection  from  the  Labrador  Sea  on  the  time  scale 
of  ~2  years,  with  little  or  no  local  convection  east  of  Greenland.  One  such  study  (Bersch 
et  al,  2007)  presents  evidence  that  there  has  been  an  intermediate  salinity  maximum  layer 
continuously  present  in  the  Irminger  Sea  since  1991,  implying  that  deep  overturning  could 
not  have  occurred  during  this  time  period.  In  the  present  section  we  examine  this  issue 
using  the  collection  of  AR7E  sections.  To  further  elucidate  the  issue  of  local  versus  remote 
ventilation,  we  use  a  volumetric  T-S  approach  applied  by  Yashayaev  (2007)  to  the  Labrador 
Basin.  We  find  strong  evidence  that  deep  convection  occurred  in  the  Irminger  Basin  dur¬ 
ing  the  early  1990s.  Furthermore,  we  identify  a  separate  signal  of  remotely-formed  LSW 
advected  into  the  Irminger  Sea  from  the  west. 

4.7.1  Intermediate  salinity  maximum  layer 

One  of  the  pieces  of  evidence  used  to  discount  the  notion  of  local  convection  in  the  Irminger 
Sea  stems  from  the  recent  analysis  of  Bersch  et  al  (2007).  They  used  time  series  obtained 
from  a  subset  of  the  AR7E  hydrographic  sections  averaged  between  34-4 1°W  (i.e.  the  full 
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width  of  the  deep  Irminger  Basin,  Fig.  4-3)  to  investigate  changes  in  the  characteristics 
and  distribution  of  LSW  in  the  Irminger  Sea.  Computing  the  temporal  changes  of  salinity 
relative  to  1991  (A591),  they  found  a  layer  of  positive  A5'9i  at  a  density  of  approximately 
Oi  =  32.38  kgjn?  (ai.5  =  34.655  kgjn? ,  between  600  and  1200  dV)  that  was  continuously 
present  since  1991  (Fig.  4-22,  right).  This  layer  coincides  with  the  Icelandic  Slope  Water 
intermediate  salinity  maximum  that  separated  the  lower  and  upper  LSW  water  masses  of 
the  2001-05  pentad  (Fig.  4-12d,  see  also  Holliday  et  al.,  2009).  Bersch  et  al.  (2007)  inter¬ 
preted  this  as  evidence  that  there  was  no  significant  formation  of  the  densest  class  of  LSW 
in  the  Irminger  Sea,  even  during  the  high-NAO  period  of  the  early  1990s.  They  argued  that, 
had  there  been  convection,  fresh  surface  water  would  have  penetrated  the  Oi  =  32.38  kg/rn^ 
isopycnal  causing  an  interruption  of  the  positive  A591  layer. 


Figure  4-22:  Temporal  changes  (left)  of  salinity  (shaded)  and  density  (Oi,  contours)  and  (right) 
salinity  relative  to  1991  (positive  values  are  shaded)  on  density  surfaces  (01.5).  Both  figures  are 
based  on  data  averaged  between  33-41°W,  and  have  been  copied  from  Bersch  et  al.  (2007). 

It  is  important  to  note  that  Bersch  et  al.  (2007)  averaged  properties  over  the  entire  deep 
Irminger  Basin,  from  the  base  of  the  Greenland  continental  slope  to  the  western  flank  of 
the  Reykjanes  Ridge  (Fig  4-3).  However,  if  deep  convection  occurs  in  the  Irminger  Sea, 
it  is  likely  to  take  place  primarily  within  the  southwestern  portion  of  the  basin  (Pickart 
et  al,  2003a,b;  Centurioni  and  Gould,  2004).  This  region  is  under  direct  influence  of  the 
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cyclonic  Irminger  Gyre  and  the  Greenland  tip  jet,  and  is  the  part  of  the  Irminger  Sea  where 
the  oceanic  and  atmospheric  conditions  are  most  favorable  for  deep  overturning.  With  this 
in  mind,  we  re-did  the  calculation  of  Bersch  et  al.  (2007)  using  hydrographic  data  confined 
to  the  location  of  the  IG  along  AR7E  (150-350  km,  see  Fig  3-lOc;  the  results  are  not  very 
sensitive  to  the  exact  limits). 

A  different  result  from  Bersch  et  al.  (2007)  is  obtained  when  the  data  are  confined 
geographically  as  such.  Figure  4-23a  shows  the  temporal  evolution  of  salinity  and  potential 
density.  In  agreement  with  Bersch  et  al.  (2007),  the  peak  of  the  FSW  dome  occurred  in 
1995.  By  2001  it  had  collapsed  (see  also  Fig.  4-12d).  However,  when  considering  only  the 
western  portion  of  the  Irminger  Basin  within  the  IG,  the  isopycnal  doming  was  much  more 
pronounced  (compare  Figures  4-22,  left,  and  4-23a).  In  particular,  the  O1.5  =  34.63  kg/rn^ 
contour  that,  according  to  Bersch  et  al.  (2007),  delimits  the  density  of  locally  formed  water 
masses  reached  a  minimum  pressure  of  approximately  230  db  in  1995.  Keep  in  mind  that 
the  data  used  here  were  collected  in  summer.  It  is  highly  unlikely  that  the  mixed  layer  depth 
in  this  region  did  not  exceed  230  db  during  the  winter  of  1994-95  (during  the  high-NAO 
period).  This  suggests  strongly  that  the  IFSW  layer  was  locally  ventilated.  Mixed  layers 
deeper  than  400  db  were  routinely  observed  by  profiling  floats^  and  moorings  during  the 
winters  of  the  early  2000s,  which,  unlike  the  early  1990s,  were  characterized  by  a  low  NAO 
index  and  a  well- stratified  water  column  (Centurioni  and  Gould,  2004;  Vage  et  al,  2008, 
2009a;  Fouarn  et  al,  2009).  Furthermore,  our  calculation  of  A591  using  hydrographic  data 
from  within  the  IG  does  indeed  show  an  interruption  of  the  positive  A591  layer  between 
1994  and  1996,  which  again  suggests  that  IFSW  was  formed  during  this  period.  This 
calculation  reinforces  the  arguments  made  to  date  that,  under  strong  atmospheric  forcing, 
deep  convection  can  occur  in  the  Irminger  Sea.  It  also  illustrates  some  of  the  subtleties 
involved  in  investigating  convection  in  that  basin. 

4.7.2  Labrador  Sea  Water  layer  thickness 

Encouraged  by  the  above  results  showing  that  deep  convection  likely  occurred  in  the  west¬ 
ern  Irminger  Sea  in  the  early  1990s,  but  realizing  that  FSW  was  also  likely  adverted  out  of 

*Argo  floats  were  not  present  in  the  early  1990s. 
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Figure  4-23:  As  in  Figure  4-22,  but  based  on  data  from  within  the  IG.  The  dashed  lines  indicate 
when  measurements  were  made.  The  ai  .5  =  34.63  contour  is  highlighted  in  red. 


the  Labrador  Sea  to  the  east,  we  studied  in  greater  detail  the  evolution  of  the  LSW  inside 
and  outside  of  the  Irminger  Gyre.  The  LSW  core  was  identified  following  the  volumetric 
density  layer  approach  of  Yashayaev  (2007).  The  method  involves  calculation  of  the  thick¬ 
ness  of  a  continuum  of  potential  density  (02)  layers.  A  mean  vertical  density  profile  for 
each  section  within  the  chosen  lateral  boundaries  was  first  determined.  The  vertical  extent 
(thickness)  of  every  layer  of  width  Aa2  =  0.01  kg/rn^  was  then  computed.  Individual  layers 
were  separated  by  0.002  kg/n? .  The  LSW  core  was  then  easily  identified  by  the  density 
layer  of  largest  thickness.  An  advantage  of  this  method  is  its  ability  to  take  into  account 
changes  in  water  mass  properties  since  it  does  not  rely  on  time-invariant  hydrographic 
bounds  to  identify  the  water  mass. 

The  spreading  of  LSW  from  the  Labrador  Sea  was  investigated  by  Pickart  et  al.  (2003b) 
using  the  advective-diffusive  model  of  Straneo  et  al.  (2003),  which  was  based  on  a  realistic 
flow  field  determined  from  float  measurements  containing  the  interior  recirculations  in  the 
Labrador  and  Irminger  seas  (Lavender  et  al,  2000).  In  agreement  with  recent  observa- 
tionally  based  estimates  (e.g.  Yashayaev,  2007;  Falina  et  al.,  2007),  the  time  for  LSW  to 
travel  from  the  Labrador  Basin  to  the  central  Irminger  Basin  was  approximately  2  years. 
Once  inside  the  Irminger  Sea,  the  model  tracer  was  observed  to  wrap  around  the  Irminger 
Gyre  (Straneo  et  al.,  2003),  taking  another  year  to  penetrate  to  the  gyre’s  center  (Pickart 
et  al.,  2003b).  The  LSW  formed  in  the  Labrador  Basin  is  in  general  colder,  fresher,  and 
denser  than  the  LSW  formed  locally  in  the  Irminger  Sea  (Pickart  et  al,  2003b).  Hence 
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we  expect  that  the  arrival  of  the  remote  LSW  will  present  an  unmistakable  signal  in  the 
Irminger  Sea  LSW  timeseries  derived  from  the  volumetric  approach.  Furthermore,  the 
strongest  remotely-formed  LSW  signal  should  occur  outside  of  the  gyre  (i.e.  to  the  east  of 
the  northward  limb  of  the  gyre). 

As  a  first  calculation  we  applied  the  volumetric  method  to  the  center  of  the  Irminger 
Gyre,  between  200  and  300  km  along  AR7E.  These  narrow  bounds  (on  average  the  gyre 
was  located  between  150  and  350  km)  were  chosen  to  avoid  contamination  from  LSW 
outside  of  the  gyre,  but  the  outcome  is  not  very  sensitive  to  this  particular  choice.  The 
resulting  average  layer  thickness  versus  time  and  density  is  shown  in  Figure  4-24. 
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Figure  4-24:  Average  thickness  (in  m)  of  potential  density  layers  within  the  Irminger  Gyre  (100- 
200  km  along  AR7E)  in  relation  to  potential  density  {02,kg/m^)  of  the  given  layer  and  time  (years). 
The  vertical  extent  of  each  overlapping  layer  was  0.01  kg/rn^,  and  individual  layers  were  separated 
by  0.002  kg/rn^.  The  red  line  follows  the  crest  of  the  thickness  ridge  and  hence  identifies  the  ILSW 
core  in  each  section.  The  dashed  lines  indicate  when  individual  sections  were  made. 


The  most  notable  feature  in  Figure  4-24  is  the  LSW  dome  that  reached  its  greatest  vol¬ 
ume  in  1995  (in  good  agreement  with  Figure  4-23a  and  Bersch  et  al,  2007).  As  previously 
mentioned,  this  agrees  well  with  the  modeled  mixed  layer  depth  in  the  winter  of  1994-95 
exceeding  1700  m  (Vage  et  al,  2008).  Indeed,  the  bottom  of  the  thickest  1995  02  level  was 
found  at  a  depth  of  1730  m  in  summer  1995  (Fig.  4-25b).  In  contrast,  the  densest  water 
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in  the  ILSW  class  was  formed  in  the  Labrador  Sea  in  the  winter  of  1993-94  (Yashayaev, 
2007),  although  ILSW  formation  did  not  shut  down  until  after  the  following  winter  (e.g 
Avsic  et  ah,  2006).  Had  the  LSW  observed  inside  the  Irminger  Gyre  been  ventilated  re¬ 
motely  in  the  Labrador  Sea,  one  would  expect  to  see  the  thickest  LSW  layer  at  the  center 
of  the  gyre  at  the  earliest  three  years  later,  in  1997.  At  this  time,  however,  the  Irminger 
Gyre  ILSW  dome  was  already  collapsing.  Timeseries  of  ILSW  core  properties  within  the 
Irminger  Gyre  are  shown  in  Fig.  4-25.  One  sees  a  single  extremum  in  the  layer  thickness, 
salinity,  and  temperature.  (We  note  that  even  though  the  layer  thickness  started  declining 
immediately  in  1996,  the  temperature  and  salinity  remained  low  for  one  more  year,  perhaps 
due  to  some  amount  of  remotely-formed  LSW  penetrating  the  gyre.) 


Figure  4-25:  Timeseries  of  core  ILSW  properties,  a)  Layer  thickness  (m).  b)  Upper  and  lower 
bounds  of  the  layer  in  a),  c)  Salinity,  d)  Potential  temperature  (°C).  One  standard  deviation  is 
indicated  by  the  error  bars. 


Note  also  the  arrival  of  the  shallower  LSW  mode  around  2000.  While  not  as  prominent 
as  in  the  Labrador  Sea  (see  Figure  6a  in  Yashayaev,  2007),  some  uLSW  evidently  formed 
also  within  the  Irminger  Gyre.  The  thickness  of  the  uLSW  layer  appears  to  have  increased 
substantially  between  2004  and  2005,  perhaps  due  to  arrival  of  uLSW  from  the  Labrador 
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Sea. 

Now  we  contrast  the  results  from  the  volumetric  density  layer  calculation  within  the 
Irminger  Gyre  to  that  outside  of  the  gyre.  Figure  4-26  corresponds  to  Figure  4-24,  but  is 
obtained  using  data  from  the  eastern  part  of  the  Irminger  Basin,  between  400  and  500  km 
along  AR7E.  The  bounds  were  chosen  such  that  the  domain  was  well  outside  the  Irminger 
Gyre,  while  still  to  the  west  of  the  Irminger  Current.  The  difference  between  the  two  figures 
is  striking.  The  ILSW  dome  is  much  smaller  in  magnitude  (the  maximum  ILSW  thickness 
was  about  500  m  compared  to  more  than  800  m  inside  the  gyre),  but  now  there  is  a  second 
peak  in  1997.  This  is  seen  more  clearly  in  the  timeseries  of  layer  thickness,  salinity,  and 
temperature  outside  of  the  gyre  (Fig.  4-27).  In  contrast  to  the  center  of  the  gyre,  there  are 
now  two  extrema  in  these  properties.  The  timing  of  the  second  peak  (1997)  is  consistent 
with  a  2-3  year  time  lag  from  the  Labrador  Sea.  Furthermore,  the  second  peak  is  colder, 
fresher,  and  deeper,  consistent  with  the  more  extreme  properties  of  the  LSW  formed  in 
the  Labrador  Basin.  Linally,  the  mean  absolute  velocity  section  for  the  pentad  of  the  early 
1990s  (Lig  4- 13c)  shows  northeastward  flow  in  the  region  outside  of  the  gyre  (between  400 
and  500  km),  consistent  with  a  remote  pathway  from  the  Labrador  Sea.  We  note  that  the 
earlier  LSW  peak  implies  some  degree  of  local  convection  outside  of  the  gyre,  which  was 
observed  in  winter  2007-08  (Lemke  de  Jong,  personal  communication,  2008). 

To  summarize,  we  have  demonstrated  that  inside  of  the  Irminger  Gyre  there  was  no 
mid-depth  salinity  maximum  layer  that  indicated  a  limit  to  the  depth  of  local  convection 
in  the  early  1990s.  Instead  there  was  a  single  peak  of  ILSW  in  1995,  too  early  to  be  a 
remote  signal  from  the  maximum  convection  year  in  the  Labrador  Sea  (1994).  However, 
outside  of  the  gyre  there  were  two  maxima.  The  timing  and  properties  of  the  latter  peak  are 
indeed  only  reconcilable  with  a  remote  Labrador  Sea  source,  and  consistent  as  well  with 
the  velocity  field. 

4.8  Discussion 

We  found  above  that  the  Irminger  Gyre  circulation  had  intensified  between  the  early  1990s 
and  the  early  2000s,  with  the  transport  increasing  at  a  rate  of  4.3  Sv  per  decade.  This 
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Figure  4-26:  As  Figure  4-24,  but  for  the  region  east  of  the  Irminger  Gyre  between  400  and  500  km 
along  AR7E. 
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Figure  4-27:  As  Figure  4-25,  but  for  the  region  east  of  the  Irminger  Gyre  between  400  and  500 
km  along  AR7E.  The  timeseries  were  cut  before  2001  since  reliable  detection  of  the  ILSW  was  no 
longer  feasible  (the  noisy  2000s  signal  is  evident  in  Figure  4-26). 


114 


was  unexpected  given  the  general  decline  in  atmospheric  forcing  that  occurred  through  the 
period  of  observations.  Inspection  of  the  AVISO  dynamic  ocean  topography  field  revealed 
that  the  general  increase  in  SSH  across  the  entire  subpolar  gyre  was  not  uniform  within 
the  Irminger  Sea,  but  that  the  eastern  half  of  the  basin  rebounded  to  a  greater  degree  than 
the  western,  leading  to  a  steeper  gradient  across  the  IG  and  hence  the  stronger  circulation 
(Fig.  4-16).  Analysis  of  the  AR7E  hydrographic  data  indicate  that  part  of  the  SSH  change 
was  due  to  steric  height  differences.  In  particular,  from  the  early  1990s  to  the  early  2000s 
the  steric  height  in  the  eastern  part  of  the  Irminger  Basin  increased  more  rapidly  than  in 
the  western  part,  consistent  with  the  evolution  of  absolute  SSH.  This  suggests  that  the 
restratification  process  leading  to  the  collapse  of  the  ILSW  dome  was  more  pronounced  on 
the  eastern  side  of  the  basin,  which  is  not  surprising  due  to  the  presence  of  the  buoyant 
Irminger  Current  next  to  the  Reykjanes  Ridge. 

As  ILSW  formation  ceased  and  the  subpolar  gyre  contracted  through  the  1990s  (e.g. 
Hakkinen  and  Rhines,  2004),  increased  flow  of  warm,  saline  subtropical  water  into  the 
eastern  part  of  the  North  Atlantic  ensued  (Hatun  et  ah,  2005).  This  in  turn  affected  the 
hydrographic  properties  of  the  IC.  In  order  to  influence  the  interior  Irminger  Sea,  the 
Atlantic-origin  water  must  be  fluxed  out  of  the  boundary  current.  Lateral  exchange  be¬ 
tween  the  boundary  current  and  the  interior,  typically  in  the  form  of  mesoscale  eddies,  has 
been  extensively  investigated  in  the  eastern  Labrador  Sea  (e.g.  Lilly  et  ai,  2003;  Straneo, 
2006).  A  similar  process  likely  takes  place  in  the  eastern  Irminger  Sea  as  well.  A  second 
mechanism  leading  to  restratification  is  transport  by  baroclinic  Rossby  waves.  In  the  sub¬ 
polar  North  Atlantic  the  internal  deformation  radius  is  on  the  order  of  10  km  (Chelton  et  al, 
1998).  Accordingly,  baroclinic  Rossby  waves  would  be  able  to  cover  the  distance  between 
the  IC  and  the  IG  in  less  than  10  years.  Though  we  did  not  see  clear  evidence  of  such 
a  westward-propagating  signal  in  the  hydrographic  data,  it  could  be  masked  by  synoptic 
variability  or  the  competing  effect  of  a  cold,  fresh  influx  along  the  spreading  pathway  of 
LSW  from  the  Labrador  Sea.  If  westward  propagation  was  the  primary  mechanism  of  the 
IG  intensification,  the  assumption  of  mass  balance  between  the  northward  and  southward 
branches  of  the  gyre  may  be  called  into  question.  In  this  case  water  masses  to  the  southeast 
of  the  originally  defined  northward  limb  of  the  gyre  may  flow  northward  alongside  the  gyre 
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but  not  recirculate  within  the  gyre. 

Even  though  intensification  of  the  IG  circulation  through  a  period  of  diminishing  wind 
stress  curl  appears  paradoxical,  there  is  not  necessarily  a  contradiction  between  our  find¬ 
ings  and  the  model  results  of  Spall  and  Pickart  (2003).  Their  model  gyre  was  exposed  to 
wind  forcing  during  winter  only,  yet  the  circulation  was  maintained  with  minimal  seasonal 
variability.  An  important  aspect  in  this  was  the  slow  propagation  speed  of  the  baroclinic 
Rossby  waves,  which  results  in  a  long  response  time  of  the  Irminger  Gyre.  As  such,  one 
would  expect  the  gyre  to  start  spinning  down  as  it  slowly  adjusts  to  the  present  reduced 
level  of  forcing.  Clearly  this  needs  further  investigation. 

The  disappearance  of  the  potential  vorticity  (PV)  anomaly  between  the  early  1990s  and 
the  early  2000s  (Fig.  4-13)  provides  an  opportunity  to  estimate  the  diffusion  coefficient 
directly  from  observations  using  the  diffusion  equation 

Dx  9 

—  =  KV\  (4.1) 

Dt 

where  K  is  the  diffusion  coefficient  and  x  is  the  tracer  quantity  of  interest,  in  this  case 
PV.  The  change  in  PV  at  the  LSW  level  within  the  IG  over  the  10-year  period  and  the 
Laplacian  of  the  PV  field  between  the  IG  and  the  IC  were  approximated  from  Figure  4- 
13,  resulting  in  a  diffusion  coefficient  of  about  160  /s.  For  comparison,  the  diffusion 

coefficient  was  also  evaluated  using  the  parametric  representation  of  Visbeck  et  al.  (1997). 
Guided  by  baroclinic  instability  theory,  Visbeck  et  al.  (1997)  proposed  a  transfer  coefficient 
(corresponding  to  the  diffusion  coefficient)  of  the  form 

K  =  a-^f  (4.2) 

yRi 

where  a  is  a  proportionality  constant  (empirically  found  to  be  0.015),  /  is  the  Coriolis 
parameter,  Ri  =  N'^ /{dvjdz)^  is  the  Richardson  number  of  the  large-scale  flow,  and  I  is 
a  length  scale  measuring  the  width  of  the  baroclinic  zone.  The  mixing  length  scale  was 
estimated  from  the  Rhines  scale  {  \JU /|3,  Rhines,  1975)  using  a  characteristic  current  speed 
of  0.1  m/s,  resulting  in  approximately  90  km  in  the  subpolar  North  Atlantic.  The  quantity 
f/V^  is  the  growth  rate  of  an  Eady  instability  (Eady,  1949),  and  the  physical  basis  for 
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the  transfer  coefficient  is  the  product  of  an  eddy  velocity  fl/y/Wi  and  a  length  scale  1. 
As  such,  the  mean  Visbeck  et  al.  (1997)  transfer  coefficient  in  the  IG  was  180  ni^/s,  and 
approximately  one  third  of  this  value  if  the  LSW  layer  is  considered  separately.  This  result 
is  in  surprisingly  good  agreement  with  the  coefficient  estimated  from  the  diffusion  equation 
(Eq.  4.1),  but  below  the  range  of  300-2000  m^/v  calculated  by  Visbeck  et  al.  (1997).  They 
found  the  minimum  value  in  a  convective  chimney,  under  conditions  very  similar  to  the 
weakly  stratified  LSW  layer  considered  here.  The  good  agreement  between  the  coefficients 
stemming  from  the  diffusion  equation  and  the  parametric  representation  may  suggest  that 
the  form  proposed  by  Visbeck  et  al.  (1997)  is  a  reasonable  representation  of  eddy  transfer, 
assuming  that  eddy  diffusion  was  the  primary  mechanism  leading  to  the  collapse  of  the 
LSW  dome.  It  is,  however,  unlikely  that  advection  did  not  play  a  role  as  well. 

A  structural  change  in  the  deep  Irminger  Current  was  observed  between  the  early  1990s 
and  the  early  2000s  (Lig.  4-13).  In  the  first  pentad  there  was  primarily  a  single  deep  ex¬ 
tension  of  the  IC  located  near  the  crest  of  the  Reykjanes  Ridge,  while  in  the  latter  pentad  a 
second  deep  extension  of  the  current  emerged  farther  down  the  ridge  near  the  2500  m  iso¬ 
bath.  This  latter  feature  carried  the  bulk  of  the  deep  IC  transport  in  the  early  2000s,  and  was 
chiefly  responsible  for  the  overall  increase  in  top  to  bottom  IC  transport  noted  above.  The 
hydrographic  properties  of  the  two  extensions  differ  (this  is  evident  in  Ligure  4-12).  While 
still  of  Atlantic  origin,  the  water  transported  by  the  western  deep  extension  is  colder  and 
fresher  than  that  in  the  eastern,  which  suggests  that  the  two  deep  extensions  have  different 
origins. 

Such  a  notion  is  consistent  with  synoptic  hydrographic  surveys  of  the  eastern  North 
Atlantic  (e.g.  Read,  2001;  Pollard  et  al.,  2004)  and  subsurface  float  data  (Bower  et  al., 
2002).  These  studies  have  revealed  that  there  are  two  distinct  pathways  from  the  North 
Atlantic  Current  to  the  Irminger  Current.  The  outer  pathway  follows  the  boundary  of  the 
subpolar  gyre  around  the  rim  of  the  Iceland  Basin,  while  the  inner  pathway  bifurcates  from 
the  North  Atlantic  Current  near  30°W  and  flows  northward,  bypassing  the  Iceland  basin. 
Near  the  southeastern  corner  of  the  Irminger  Basin  the  two  pathways  merge  to  form  the 
Irminger  Current  and  proceed  jointly  northwards  along  the  western  flank  of  the  Reykjanes 
Ridge.  Using  hydrographic  data,  Pollard  et  al.  (2004)  demonstrated  that  these  pathways 
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shifted  between  1991  and  1996,  apparently  due  to  the  contraction  of  the  subpolar  gyre.  It  is 
not  unreasonable  to  think  that  the  relative  strengths  of  these  two  pathways  may  change  in 
time  as  well,  and  that  the  increasing  flow  observed  in  the  western  IC  extension  could  be  a 
reflection  of  the  large-scale  changes  that  took  place  throughout  the  subpolar  gyre  from  the 
early  1990s  to  the  early  2000s. 

4.9  Summary  and  conclusions 

Significant  changes  took  place  in  the  subpolar  North  Atlantic  between  the  early  1990s  and 
early  2000s.  Associated  with  a  significant  weakening  in  the  atmospheric  forcing,  the  for¬ 
mation  of  lower  LSW  in  the  Labrador  Sea  all  but  ceased,  and  the  subpolar  gyre  circulation 
slowed  and  contracted. 

In  this  study  we  analyzed  a  collection  of  32  hydrographic  transects  across  the  southern 
Irminger  Sea  from  the  period  1991  to  2007.  The  dramatic  changes  that  took  place  in  the 
broad  subpolar  region  are  reflected  in  the  Irminger  Sea  hydrographic  data.  Most  notable  is 
the  collapse  of  the  lower  LSW  dome  that  reached  peak  volume  in  1995.  An  EOF  analysis 
identified  these  changes  as  the  primary  mode  of  variability,  and  demonstrated  a  high  degree 
of  coherence  between  the  subpolar  gyre  index  of  Hakkinen  and  Rhines  (2004)  calculated 
from  North  Atlantic  SSH  levels  and  the  evolution  of  the  local  Irminger  Sea  hydrography. 

From  absolute  geostrophic  velocity  sections  referenced  using  satellite  altimetry,  we 
quantified  for  the  first  time  some  of  the  characteristics  of  the  Irminger  Gyre,  a  closed, 
cyclonic  recirculation  in  the  western  Irminger  Sea.  The  mean  transport  of  the  gyre  was  6.8 
±  1.8  Sv,  and  over  the  17-year  study  period  the  transport  of  the  gyre  increased  at  a  rate  of 
4.3  Sv  per  decade.  This  increase  was  associated  with  a  steeper  sea  level  gradient,  which  is 
partially  due  to  changes  in  steric  height.  In  particular,  buoyant  water,  presumably  fluxed 
from  the  Irminger  Current,  enhanced  the  restratification  process  in  the  eastern  Irminger 
Basin.  This  is  consistent  with  the  increased  presence  of  subtropical-origin  water  in  the 
Irminger  Current  during  the  second  half  of  the  study  period,  related  to  the  contraction  of 
the  subpolar  gyre. 

Consistent  with  the  notion  of  a  declining  subpolar  gyre,  the  upper  layer  transport  of 
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the  Irminger  Current  decreased  over  the  study  period.  However,  the  full-depth  transport  of 
the  current  increased  due  to  the  emergence  of  a  second  deep  extension  of  the  IC.  Relative 
changes  between  two  separate  pathways  that  contribute  to  the  IC  may  have  affected  this 
development.  The  mean  northward  full-depth  IC  transport  was  estimated  to  20  ±  3  Sv. 
On  the  western  side  of  the  basin  the  mean  net  southward  full-depth  Western  Boundary 
Current  transport  was  -32  ±  4  Sv  (this  is  the  equatorward  throughput,  as  the  southward 
limb  of  the  IG  has  been  subtracted  out).  For  the  DWBC  (denser  than  Oe  >  27.8  kg/m?) 
we  determined  a  mean  net  transport  of  -9.6  ±1.3  Sv,  in  good  agreement  with  previous 
estimates.  Consistent  with  other  studies  in  both  the  Irminger  and  Labrador  seas  we  found 
that  the  DWBC  flow  increased  through  the  study  period. 

Recognizing  the  importance  of  sub-basin  scale  variation  in  the  atmosphere  and  ocean 
in  regard  to  convection  in  the  Irminger  Sea,  we  examined  the  hydrographic  data  both  in¬ 
side  and  outside  of  the  Irminger  Gyre,  respectively.  We  demonstrated  that  the  intermediate 
salinity  maximum  purportedly  placing  a  limit  on  the  depth  of  convection  in  the  Irminger 
Sea  (Bersch  et  al,  2007)  was  not  continuously  present  within  the  IG.  Furthermore,  through 
a  volumetric  analysis  of  the  ILSW  layer  within  the  Irminger  Gyre,  we  found  a  single  max¬ 
imum  in  layer  thickness  in  1995,  too  early  to  be  a  remote  signal  from  the  Labrador  Sea.  In 
contrast,  outside  of  the  gyre  there  were  two  ILSW  layer  thickness  peaks.  The  first  (in  1994) 
was  likely  caused  by  local  convection,  while  the  second  (in  1997)  is  only  reconcilable  in 
timing  and  properties  with  a  remote  Labrador  Sea  source.  We  believe  that  our  results  thus 
clarify  the  relative  roles  of  local  versus  remote  convection  in  contributing  to  the  presence 
of  LSW  in  the  Irminger  Sea. 
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Chapter  5 


Conclusions  and  Discussion 


In  this  thesis  deep  convection  and  Labrador  Sea  Water  (LSW)  formation  in  the  subpo¬ 
lar  North  Atlantic  have  been  investigated,  with  special  emphasis  on  the  Irminger  Sea  and 
the  small-scale  atmospheric  and  oceanic  processes  influencing  local  overturning.  Some 
of  the  main  results  are:  (i)  small-scale  atmospheric  (the  westerly  Greenland  tip  jet)  and 
oceanic  (the  Irminger  Gyre)  features  are  of  fundamental  importance  to  deep  convection  in 
the  Irminger  Sea,  (ii)  deep  mixed  layers  in  the  Irminger  Sea  were  for  the  first  time  docu¬ 
mented  in  winter  2007-08,  and  the  reasons  why  the  overturning  took  place  were  elucidated, 
and  (iii)  the  competing  signals  of  local  versus  remote  formation  of  LSW  were  clarified  and 
distinguished,  suggesting  that  LSW  was  formed  in  the  Irminger  Sea  in  the  early  1990s.  ^ 

5.1  Ramifications  of  LSW  formation  in  the  Irminger  Sea 

If  LSW  is  indeed  formed  outside  of  the  Labrador  Sea,  this  alters  our  current  understand¬ 
ing  of  the  convective  process  in  the  subpolar  North  Atlantic.  Presently,  there  exists  in  the 
literature  a  wide  range  of  LSW  formation  rates  (e.g.  Pickart  and  Spall,  2007;  Haine  et  al, 
2008),  due  in  part  to  different  methods  of  estimation  and  definitions  of  LSW,  as  well  as 
the  variability  of  convective  activity.  This  has  made  assessment  of  LSW  production  prob¬ 
lematic.  An  additional  source  of  LSW  further  compounds  this  difficulty,  making  it  harder 
to  quantify  other  aspects  of  LSW  formation  such  as  spreading  rates  and  ventilation  times. 

^Please  refer  to  Chapters  2-4  for  more  detailed  conclusions  and  discussions. 
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The  notion  of  LSW  formation  in  the  Irminger  Sea  entails  a  new  perspective  in  which  to 
interpret  both  observations  and  model  results,  including  a  reassessment  of  the  influence  of 
LSW  on  NADW  formation  and  the  Meridional  Overturning  Circulation  (MOC). 

Another  potential  ramification  is  that  convection  in  the  Irminger  Sea  could  have  a  more 
direct  impact  on  the  weather  and  climate  in  northwestern  Europe  than  convection  in  the 
Labrador  Sea.  Given  its  closer  proximity  to  the  North  Atlantic  storm  track  and  the  Icelandic 
Low,  release  of  heat  through  convection  in  the  southwestern  Irminger  Sea  may  to  a  greater 
extent  energize  the  cyclones  that  subsequently  track  towards  Europe.  On  the  other  hand, 
proper  resolution  of  tip  jets,  which  are  responsible  for  the  greatest  heat  fluxes  from  the 
Irminger  Sea,  does  not  seem  to  have  a  large  impact  on  European  weather  in  a  numerical 
weather  prediction  model  (Jung  and  Rhines,  2007). 

5.2  Comparison  with  hindcast  simulations 

In  this  section,  difficulties  associated  with  modeling  the  subpolar  North  Atlantic  are  dis¬ 
cussed.  The  results  of  recent  hindcast  simulations  compared  to  observations  follows,  with 
a  special  emphasis  on  convection  in  the  Irminger  Sea.  Einally  some  suggestions  to  improve 
the  realism  of  the  simulations  are  provided. 

High-resolution  general  circulation  model  (GCM)  hindcasts  face  difficult  challenges 
in  the  subpolar  North  Atlantic.  In  particular,  small  hydrographic  contrasts  between  wa¬ 
ter  masses  are  dynamically  important,  but  difficult  to  maintain  (e.g.  Treguier  et  al,  2005; 
Haine  et  al,  2008).  Lew  models  realistically  reproduce  the  observed  hydrographic  prop¬ 
erties  (e.g.  Treguier  et  al,  2005;  de  Jong  et  al,  2009),  and  it  was  recently  found  that  the 
GCM  used  by  Yaeger  and  Jochum  (2009)  was  unable  to  accurately  simulate  characteris¬ 
tics  of  both  water  mass  properties  and  gyre  circulation  simultaneously,  which  may  have  a 
bearing  also  on  other  simulations.  Indeed,  most  of  the  recent  high-resolution  GCM  runs  in¬ 
vestigated  by  Treguier  et  al  (2005)  did  produce  a  realistic  circulation,  but  were  not  close  to 
the  observed  temperature  and  salinity  fields.  In  the  subpolar  North  Atlantic,  a  realistic  sim¬ 
ulation  of  the  hydrographic  properties  takes  on  greater  importance  through  the  formation 
of  LSW  and  its  impact  on  climate  and  the  MOC. 
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Excessive  convection  is  a  common  problem  to  many  models,  and  tends  to  result  from 
a  lack  of  restratification  (Treguier  et  al,  2005;  Maine  et  al,  2008;  de  Jong  et  al,  2009). 
Direct  simulation  of  eddy  fluxes  through  further  increase  in  resolution  is  likely  the  only 
way  to  overcome  this  issue.  The  Nordic  seas  overflows,  in  particular  their  magnitude  and 
the  amount  of  entrainment  and  mixing  associated  with  their  descent  from  the  Greenland- 
leeland-Scotland  sills,  are  also  currently  difficult  to  realistically  reproduce.  Very  high  hor¬ 
izontal  and  vertical  resolution  is  required,  in  particular  in  the  region  where  the  overflows 
descend  and  are  modified  by  ambient  waters.  However,  this  is  not  yet  possible  with  the 
currently  available  computational  power  (Maine  et  al,  2008).  This  undoubtedly  has  an  ef¬ 
fect  on  the  realism  of  the  hindcasts,  and  may  alter  the  sensitivity  of  the  MOC  to  the  rate  of 
LSW  formation  in  the  subpolar  gyre  common  to  most  simulations. 

A  set  of  recent  hindcast  simulations  (e.g.  Hatun  et  al,  2005;  Treguier  et  al,  2005; 
Boning  et  al,  2006;  Deshayes  and  Frankignoul,  2008)  have  reproduced  the  overall  decline 
of  the  subpolar  gyre  eirculation  documented  by  Hakkinen  and  Rhines  (2004),  as  well  as 
the  reduced  levels  of  convective  activity  since  the  mid-1990s.  This  is  consistent  with  the 
results  of  Chapter  4,  as  is  the  increase  in  temperature  and  salinity  resulting  from  greater 
volumes  of  subtropical  water  in  the  eastern  North  Atlantic  following  the  contraction  of  the 
subpolar  gyre  (Berseh  et  al,  2007;  Sarafanov,  2009).  However,  while  these  large-scale 
features  have  been  realistically  simulated,  the  modeled  hydrographic  properties  deviate 
signifieantly  from  measured  values  (Treguier  et  al,  2005;  de  Jong  et  al,  2009).  Salinity 
drift  is  a  partieularly  difficult  problem  to  solve  (e.g.  Treguier  et  al,  2005). 

As  an  explicit  example,  consider  the  location  and  depth  of  convection  in  the  western 
subpolar  gyre.  Deep  convection  in  the  Irminger  Sea  was  reported  in  several  studies  (e.g. 
Treguier  et  al,  2005;  Deshayes  et  al,  2007).  However,  close  inspection  of  the  hindcast  an¬ 
alyzed  by  Deshayes  et  al  (2007)  reveals  that  the  convection  took  place  in  the  northeastern 
Irminger  Sea,  in  the  vicinity  of  the  boundary  current.  According  to  the  observed  potential 
vorticity  structure  and  levels  of  atmospheric  forcing  (Pickart  et  al,  2003b),  this  is  not  where 
deep  convection  takes  place  in  the  Irminger  Sea.  Indeed,  the  LSW  formation  in  the  early 
1990s  investigated  in  Chapter  4  took  place  inside  the  Irminger  Gyre  in  the  southwestern 
Irminger  Sea.  This  difference  in  location  may  seem  subtle,  but  turns  out  to  have  consid- 
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erable  consequences.  Deshayes  and  Frankignoul  (2008)  found  the  MOC  to  be  impacted 
more  strongly  by  convection  in  the  Irminger  Sea  than  in  the  Labrador  Sea.  However,  in 
their  model,  convection  in  the  Irminger  Sea  took  place  adjacent  to  the  boundary  current, 
and  hence  immediately  influenced  the  DWBC  and  thereby  the  MOC.  In  contrast,  in  the 
Labrador  Sea  convection  took  place  within  a  recirculation  gyre,  where  it  was  to  a  greater 
degree  retained  and  prevented  from  affecting  the  MOC.  These  results  highlight  the  impor¬ 
tance  of  understanding  precisely  where  and  how  deep  convection  occurs  in  the  subpolar 
gyre  -  including  the  Irminger  Sea. 

In  order  to  produce  more  realistic  hindcasts  from  the  Irminger  Sea  and  more  properly 
assess  its  importance,  increased  resolution  in  both  models  and  the  atmospheric  products 
used  to  force  the  models  is  very  likely  of  key  importance.  Not  only  will  a  fully  eddy¬ 
resolving  model  be  capable  of  more  accurately  simulating  the  depth  of  convection  through 
proper  preconditioning  and  restratification  via  eddy  fluxes,  the  small  scale  features  impor¬ 
tant  to  convection  in  the  Irminger  Sea  (the  westerly  Greenland  tip  jet  and  the  Irminger 
Gyre)  will  also  be  properly  included. 

5.3  Effects  of  an  intensifying  Irminger  Gyre  circulation 

The  intensification  of  the  Irminger  Gyre  (IG)  through  the  16-year  period  1991-2007,  doc¬ 
umented  in  Chapter  4,  is  likely  to  have  had  an  influence  on  the  convective  activity  in  the 
Irminger  Sea.  The  obvious  notion  is  that  a  stronger  gyre  would  result  in  deeper  convection. 
For  instance,  a  stronger  circulation  would  tend  to  increase  the  gyre’s  ability  to  retain  water 
convected  in  previous  winters,  and  thus  enhance  the  preconditioning  for  the  following  win¬ 
ter.  In  concert  with  greater  isopycnal  doming,  the  stronger  gyre  circulation  would  hence 
be  conducive  for  deep  convection.  However,  this  scheme  does  not  agree  with  the  observed 
trend  of  LSW  formation,  whereby  deep  convection  took  place  in  the  early  1990s  during  the 
weaker  gyre  state,  while  there  was  no  LSW  formation  during  the  subsequent  intensifying 
phase  of  the  IG.  This  discrepancy  can  to  a  great  extent  be  explained  by  the  much  stronger 
buoyancy  forcing  during  the  first  pentad  of  the  1990s  (Fig.  4-14).  Furthermore,  during  the 
deep  convection  in  winter  2007-08,  the  deepest  mixed  layers  were  observed  farther  east  in 
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the  Irminger  Basin,  outside  of  the  IG  (Chapter  3  and  Femke  de  Jong,  personal  communi¬ 
cation,  2008).  Although  the  circulation  was  stronger  that  winter,  the  degree  of  isopycnal 
doming  was  much  reduced  compared  to  the  early  1990s  (Fig.  4-13).  This  suggests  that 
the  presence  of  a  recirculation  gyre  may  not  have  been  overly  important  to  convection  in 
winter  2007-08,  in  part  due  to  the  lack  of  preconditioning.  As  cyclonic  circulation  is  one  of 
the  criteria  for  open-ocean  convection  listed  by  Marshall  and  Schott  (1999),  its  importance 
to  deep  convection  in  the  Irminger  Sea  needs  to  be  investigated  further.  The  importance  of 
small-scale  features  to  deep  convection  in  the  Irminger  Sea  may  imply  that  these  criteria 
are  less  applicable  in  that  region. 

The  intensity  of  the  Irminger  Gyre  would  also  have  an  influence  on  the  spreading  of 
LSW  from  the  Irminger  Sea.  A  stronger  gyre  would  to  a  greater  extent  trap  the  water 
parcels  inside.  As  discussed  above,  this  is  conducive  to  deep  convection,  but  at  the  same 
time  would  also  limit  the  influence  of  the  convected  water  outside  of  the  gyre.  Despite  the 
inaccuracy  of  the  location  of  deep  convection  in  the  GCM  hindcast  analyzed  by  Deshayes 
and  Frankignoul  (2008),  it  serves  to  illustrate  the  difference  between  convection  inside 
and  outside  a  recirculation  gyre.  In  the  Labrador  Sea,  deep  convection  mainly  took  place 
within  the  gyre,  which  tended  to  retain  the  convected  product.  In  contrast,  in  the  Irminger 
Sea  convection  took  place  outside  the  gyre  in  the  vicinity  of  the  boundary  current,  and 
as  a  result  LSW  from  the  Irminger  Sea  was  directly  exported  via  the  DWBC,  resulting  in 
a  strong  influence  on  the  MOC  in  their  model.  Clearly  convected  waters  in  or  near  the 
boundary  currents  are  more  easily  exported  from  the  convective  region.^ 

The  distribution  of  LSW  in  the  subpolar  North  Atlantic  has  been  mapped  by  several 
authors  (e.g.  Talley  and  McCartney,  1982;  Straneo  et  al,  2003;  Yashayaev  et  al,  2007a), 
and  is  a  result  of  convective  activity  in  both  the  Labrador  and  Irminger  basins.  In  particular, 
the  presence  of  the  recirculation  gyres  leads  to  a  two-way  exchange  of  LSW  between  the 
Labrador  and  Irminger  seas  (Lavender  et  al,  2000;  Spall  and  Pickart,  2003),  which  is  in 
contrast  to  the  Labrador  Sea-centric  view  of  LSW  spreading  solely  from  the  Labrador  Sea 
into  the  Irminger  Sea.  In  fact,  Kvaleberg  et  al.  (2008)  found  that  the  observed  distribution 

^Convection  within  the  western  boundary  current  of  the  Labrador  Sea  has  been  observed  (Pickart  et  al. 
2002)  and  modeled  (Deshayes  and  Frankignoul,  2008),  but  the  properties  of  the  resulting  water  mass  are  less 
extreme  than  for  convection  in  the  interior. 
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of  LSW  is  in  better  agreement  with  sources  both  in  the  Labrador  and  Irminger  seas  than 
with  a  single  source  in  the  Labrador  Sea.  The  effects  of  LSW  formation  in  the  subpolar  gyre 
are  communicated  to  the  outside  via  spreading  in  the  DWBC  and  along  interior  pathways, 
and  through  its  impact  on  the  MOC  and  meridional  heat  transport.  Hence  the  influence  of 
LSW  formation  is  not  merely  a  regional  phenomenon  in  the  subpolar  North  Atlantic. 

5.4  Meridional  heat  transport  and  a  changing  climate 

Meridional  heat  transport  by  the  Atlantic  MOC  is  of  significant  importance  for  the  climate 
of  the  northern  hemisphere,  in  particular  northwestern  Europe  (e.g.  Vellinga  and  Wood, 
2002).  Numerous  examples  suggest  that  changes  in  the  MOC  will  have  climatic  conse¬ 
quences.  Past  abrupt  climate  change  documented  in  paleoclimatic  records  has  been  linked 
with  changes  in  ocean  circulation  (e.g.  Heinrich,  1988;  Dansgaard  et  al,  1993;  Rahmstorf, 
2002;  McManus  et  al,  2004).  These  events  appear  to  have  been  centered  in  the  subpolar 
North  Atlantic,  but  their  effects  were  felt  globally.  Such  episodes  are  not  limited  to  the  past. 
Of  more  immediate  concern  are  model  studies  suggesting  that  the  convective  activity  in  the 
subpolar  region  and  the  MOC  are  sensitive  to  anthropogenic  climate  change  (e.g.  Gregory 
et  al,  2005;  Stouffer  et  al,  2006). 

In  this  section  the  meridional  heat  transport  and  overturning  circulation  in  the  current 
climatic  state  is  first  discussed.  Various  mechanisms  likely  to  impact  the  MOC  in  a  warm¬ 
ing  climate  are  then  subsequently  considered.  The  section  concludes  with  a  discussion  of 
the  time  interval  of  observations  likely  required  to  detect  a  trend  resulting  from  a  changing 
climate,  and  an  outlook  on  future  work  in  this  region. 

The  MOC  in  the  presently  evolving  climatic  state 

Northward  heat  transport  into  the  North  Atlantic  peaks  at  1.3  PW  (1  PW  =  10^^  W, 
e.g.  Hall  and  Bryden,  1982)  in  the  subtropics.  Decomposing  the  overturning  circulation 
into  three  separate  components;  shallow,  intermediate,  and  deep,  where  the  latter  two  are 
associated  with  LSW  and  NADW  formation,  respectively,  Talley  (2003)  found  that  as  much 
as  one  third  is  accounted  for  by  LSW  formation.  In  contrast,  Pickart  and  Spall  (2007) 
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estimated  an  intermediate  northward  heat  flux  an  order  of  magnitude  smaller  than  this  from 
repeat  hydrographic  sections  across  the  Labrador  Sea.  They  found  that  the  heat  flux  resulted 
predominantly  from  the  horizontal  rather  than  the  overturning  circulation.  In  Pickart  and 
Spall  (2007),  only  the  northern  half  of  the  Labrador  Sea  was  included.  This  is  the  region 
where  greatest  convective  activity  and  largest  heat  loss  is  expected  to  occur,  but  a  significant 
amount  also  takes  place  outside  the  Labrador  Sea,  in  particular  in  the  Irminger  Sea,  which 
was  not  accounted  for  in  their  study. 

The  generally  shallower  depths  of  convection  in  the  subpolar  gyre  after  the  mid-1990s 
are  indicative  of  a  reduced  heat  loss  to  the  atmosphere,  and  would  tend  to  imply  a  decline 
in  the  poleward  oceanic  heat  transport.  The  resulting  contraction  of  the  subpolar  gyre,  on 
the  other  hand,  led  to  an  increase  in  the  temperature  of  the  Atlantic  inflow  to  the  Nordic 
seas  (Holliday  et  al,  2008),  which  would  enhance  the  heat  loss  there.  One  may  specu¬ 
late  that  there  is  a  compensation  between  reduced  heat  loss  in  the  subpolar  North  Atlantic 
and  increased  heat  loss  in  the  Nordic  seas  and  vice  versa.  The  poleward  heat  transport 
through  the  Atlantic  may  therefore  remain  relatively  steady,  but  the  distribution  of  heat 
would  change  significantly,  and  likely  have  an  impact  on  the  climate  of  northwestern  Eu¬ 
rope.  This  line  of  reasoning  assumes  that  there  is  not  heat  accumulation,  which  in  all 
probability  is  not  correct. 

Intensification  of  the  MOC  occurred  in  many  regional  GCMs  during  the  high-NAO 
period  ending  in  the  mid-1990s  (e.g.  Boning  etal,  2006;  Deshayes  and  Frankignoul,  2008). 
However,  this  is  not  the  case  in  the  ECCO  state  estimate  at  24°N  (Wunsch  and  Heimbach, 
2006).  There  are  a  number  of  possible  reasons  for  this  discrepancy.  On  interannual  time 
scales  a  lack  of  meridional  coherence  between  the  subpolar  and  subtropical  Atlantic  MOC 
was  noted  by  Bingham  et  al.  (2007),  and  convective  activity  in  the  subpolar  North  Atlantic 
is  poorly  represented  in  the  ECCO  model  (Schott  and  Brandt,  2007;  de  Jong  et  al,  2009). 
Despite  reduced  rates  of  LSW  formation  since  the  early  1990s,  the  DWBC  intensified  (see 
Chapter  4  and  Pickart  et  al,  1999;  Dengler  et  al,  2006;  Sarafanov  et  al,  2009b).  This  has 
not  been  noted  in  regional  hindcast  simulations,  which  tend  to  underestimate  the  overflow 
transport  and  entrainment,  and  may  therefore  attribute  excessive  influence  on  the  MOC  to 
LSW  formation.  Or  it  may  simply  be  that  the  MOC  is,  in  fact,  not  very  sensitive  to  the  rate 
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of  LSW  formation  in  the  subpolar  gyre  (Pickart  and  Spall,  2007). 


Mechanisms  impacting  convective  activity  and  the  MOC 

In  the  present  warming  climate  the  convective  activity  in  the  subpolar  North  Atlantic  is 
expected  to  diminish,  which  would  likely  impact  the  MOC  (e.g.  Stouffer  et  al,  2006).  Both 
surface  warming  and  freshening,  caused  by  Arctic  and  Greenland  ice  melt  and  an  intensi¬ 
fying  hydrological  cycle,  would  contribute  to  stabilize  the  water  column.  The  distribution 
of  this  freshwater  matters  greatly.  In  fact,  in  Chapter  3  it  was  shown  that  a  certain  amount 
of  freshwater  along  the  margins  of  the  Labrador  Sea  was  beneficial  to  deep  convection  in 
the  interior  in  winter  2007-08.  The  increased  stability  of  the  water  column  led  to  enhanced 
local  formation  of  ice,  which  shielded  the  boundary  currents  from  air-sea  interaction  and 
consequently  opened  a  pathway  for  unmodified  cold,  dry  continental  air  to  the  interior  con¬ 
vective  region.  During  that  winter  the  primary  influx  of  fresh  water  was  trough  the  Davis 
Strait,  impacting  mainly  the  northern  and  western  margins  of  the  Labrador  Sea.  This  was 
not  the  case,  however,  during  the  Great  Salinity  Anomaly  (GSA,  Dickson  et  al,  1988), 
where  a  much  larger  influx  of  freshwater  capped  the  interior  Labrador  Sea  with  a  low- 
salinity  layer  that  prevented  deep  convection  from  taking  place  for  several  winters  (Lazier, 
1980).  The  GSA  entered  the  Labrador  Sea  via  the  West  Greenland  Current,  which  is  be¬ 
lieved  to  be  the  dominant  pathway  of  freshwater  into  the  interior  of  the  basin  (Schmidt  and 
Send,  2007). 

Local  factors  are  not  the  only  agents  influencing  the  convective  activity  in  the  subpolar 
North  Atlantic.  The  origin  of  the  GSA  was  the  Arctic  Ocean.  A  change  in  the  atmospheric 
circulation  was  the  most  likely  cause  of  the  sudden  release  of  fresh  water  that  had  been 
accumulating  there.  A  similar  release  of  fresh  water  presently  accumulating  in  the  Arctic 
Ocean  may  lead  to  another  GSA  (Peterson  et  al,  2006).  Its  export  out  of  the  Arctic  via  the 
East  Greenland  Current  may  have  similar  consequences  for  the  convective  activity  when  it 
reaches  the  subpolar  gyre.  In  order  to  affect  interior  convection,  a  mechanism  of  exchange 
between  the  boundary  current  and  the  interior  is  required.  Enhanced,  topographically  gen¬ 
erated  eddy  activity  in  the  West  Greenland  Current  promotes  exchange  in  the  Labrador 


128 


Basin  (e.g.  Bracco  et  al,  2008),  but  no  corresponding  mechanism  has  been  identified  for 
the  Irminger  Sea,  which  may  thus  be  less  affected  by  a  new  GSA.  Remote  influence  also 
from  the  subtropics  in  the  form  of  sea  surface  temperature  and  heat  anomalies  has  been 
suggested  (e.g.  Sutton  and  Allen,  1997).  A  concern  is  that  this  source  of  variability  is 
excluded  from  regional  GCMs. 

Detection  of  trends  signifying  climate  change 

The  changes  in  the  subpolar  gyre’s  circulation,  extent,  and  hydrographic  properties  re¬ 
ported  in  this  thesis  and  elsewhere  (e.g.  Hakkinen  and  Rhines,  2004;  Bersch  et  al,  2007; 
Sarafanov,  2009)  that  have  taken  place  since  the  early  1990s  are  not  sufficient  to  infer 
a  long-term  trend  associated  with  a  change  in  global  climate.  Rather,  these  changes  are 
more  likely  associated  with  the  atmosphere’s  NAO-modulated  variability.^  The  observed 
variability  may  be  considered  a  superposition  of  the  oceanic  response  to  atmospheric  vari¬ 
ability  and  a  long-term  trend  stemming  from  a  slowly  changing  climate.  In  this  case  more 
than  one  cycle  of  prolonged  NAO  increase  or  decline  (such  as  the  one  that  took  place  be¬ 
tween  the  mid-1960s  and  the  present,  Sarafanov,  2009)  must  be  observed,  implying  a  need 
for  timeseries  of  at  least  multi-decadal  if  not  centennial  duration  (see,  for  instance,  Wun- 
sch,  1999).  Baehr  et  al.  (2008)  made  an  attempt  to  quantify  the  duration  of  observations 
required  to  reliably  detect  a  change  in  the  MOC  at  26°N  in  the  Atlantic.  They  found  that 
continuous  observations  and  repeated  hydrographic  transects  every  5  years  would  both  re¬ 
quire  a  period  of  about  90  years  to  detect  a  weakening  of  the  MOC  with  95%  reliability. 
This  may  be  a  relevant  time  scale  to  keep  in  mind  for  detection  of  a  trend  signifying  a 
change  in  climate  in  the  subpolar  North  Atlantic  as  well. 

Even  though  the  subpolar  North  Atlantic  is  among  the  most  well-observed  parts  of  the 
World  Ocean,  there  is  still  much  that  remains  uncertain  about  this  region  and  its  influence 
on  global  climate.  Most  model  simulations  display  a  delayed  coherence  between  LSW 
formation  and  the  vigor  of  the  MOC,  while  observation-based  analyses  produce  opposing 
results  (e.g.  Talley,  2003;  Pickart  and  Spall,  2007).  Focusing  on  LSW,  estimates  of  forma- 

^  A  divergence  between  the  NAO  index  and  an  index  describing  the  subpolar  gyre  circulation  over  the  last 
decade  was  noted  by  Hatun  et  al.  (2009).  The  large  difference  between  the  two  high-NAO  winters  of  2006-07 
and  2007-08  discussed  in  Chapter  3,  in  which  deep  convection  took  place  only  during  the  latter,  is  consistent 
with  that. 
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tion  rates  and  spreading  time  scales  have  not  converged,  and  the  Irminger  Sea  source  must 
be  taken  into  account.  In  particular,  the  relative  importance  of  the  Labrador  Sea  versus  the 
Irminger  Sea  as  a  source  of  LSW  should  be  determined.  Sustained  observations,  both  from 
space  and  in  situ,  will  be  of  crucial  importance.  Further  increase  in  model  resolution  should 
improve  the  realism  of  regional  hindcasts  and  state  estimates,  in  which  model  physics  and 
observations  are  combined.  As  a  particularly  climate-sensitive  part  of  the  World  Ocean, 
the  Irminger  Sea  and  the  subpolar  North  Atlantic  continue  to  warrant  close  observation. 
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